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Abstract 

A  commonly  used  metric  for  measuring  the  reflective  properties  of  a  sample  is  Direc¬ 
tional  Hemispherical  Reflectance  (DHR).  The  Infrared  Grazing  Angle  Reflectometer 
was  specifically  designed  and  built  to  measures  DHR  at  high  angles  of  incidence.  To 
facilitate  this  measurement,  a  novel  three  dimensional  multi-face  pyroelectric  detector 
is  used.  A  hemi-esllipsoidal  mirror  is  used  to  reflect  the  light  scattered  by  the  sample 
to  the  detector.  The  radiometric  model  for  this  system  is  derived  and  simulated  to 
illustrate  that  this  detector  can  theoretically  measure  better  than  98%  of  the  light 
reflected  by  the  sample,  regardless  of  the  diffusivity  of  the  sample  measured.  This 
error  is  shown  to  be  angularly  dependent,  and  cannot  be  calibrated  out  of  the  mea¬ 
surement.  The  setup  of  this  equipment,  as  well  as  the  modifications  implemented 
at  AFIT,  is  discussed,  and  the  procedure  for  calibration  is  derived.  Finally,  several 
sample  measurements  and  the  associated  difficulties  are  discussed. 
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CHARACTERIZATION  AND  MEASUREMENTS  FROM  THE  INFRARED 
GRAZING  ANGLE  REFLECTOMETER 

I.  Introduction 

During  November  of  2011,  the  Air  Force  Research  Lab  received  a  new  piece  of 
equipment  from  Northrop  Grumann  Corporation.  The  Infrared  Grazing  Angle  Re- 
flectometer  is  a  piece  of  equipment  designed  to  analyze  samples  and  report  back  their 
Directional  Hemispheric  Reflectance.  This  material  property  is  a  widely  used  metric 
in  the  LInited  States  Air  Force  to  describe  how  much  light  is  reflected  off  a  sample 
at  a  given  angle  of  incidence.  A  sample’s  DHR  is  used  in  many  different  ways,  rang¬ 
ing  from  modeling  the  infrared  signature  of  fighter  jets  to  calculating  how  much  the 
absorbance  of  a  material.  IGAR  quickly  and  accurately  calculates  this  value  for  a 
sample  for  a  wide  range  of  incident  angles,  increasing  the  measurement  capabilities 
of  the  Air  Force  of  the  LInited  States. 

DHR  is  a  ratio  of  the  power  reflected  by  a  surface  over  the  power  incident  on 
that  surface.  Making  measurements  of  the  power  reflected  over  such  a  large  area  is  a 
difficult  problem.  To  make  this  measurement,  Northrop  Grumman  used  a  ellipsoidal 
hemispherical  mirror  to  collect  the  light,  and  focus  it  to  a  detector.  Because  of 
the  angles  involved  in  the  measurement,  a  novel  five  sided  pyroelectric  detector  was 
created  to  measure  this  reflected  light.  In  order  to  prove  that  we  can  make  accurate 
measurements  of  the  DHR  of  the  sample,  we  must  first  prove  that  we  can  properly 
measure  this  reflected  power  using  this  detector. 

An  ideal  detector  can  measure  all  of  the  light  incident  upon  it.  For  an  actual 
detector,  this  is  not  the  case.  For  most  systems,  we  collect  light  at  only  a  single 
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angle  of  incidence,  meaning  that  this  inefficiency  is  a  constant,  and  can  be  calibrated 
out  of  the  system.  However,  for  this  system,  we  collect  light  over  a  range  of  angles. 
Because  the  ability  for  the  detector  to  absorb  light  varies  as  a  function  of  angle,  we 
can’t  calibrate  this  inefficiency  out  of  the  system.  However,  as  is  shown  later,  this 
inefficiency  is  minor,  and  can  in  most  cases  be  ignored. 

In  order  to  use  the  detector  to  make  DHR  measurements,  a  robust  calibration 
system  is  needed.  This  calibration  system  is  derived  and  discussed  for  this  system.  By 
using  mirrors  with  known  reflectivities,  we  can  calibrate  all  five  sides  of  the  detector, 
and  make  accurate  DHR  measurements  of  the  sample.  Some  of  these  measurements 
are  compared  to  their  theoretical  counterparts. 

1.1  Overview 

The  most  important  piece  of  IGAR  is  the  detector  used  in  measurements.  In 
order  to  understand  how  the  detector  makes  measurements,  the  theory  behind  the 
detector  needs  to  be  understood,  which  motivates  Chapter  III,  Detector  Theory.  This 
section  makes  two  critical  assumptions  in  the  operation  of  the  detector,  which  are  not 
true  in  practice.  Chapter  IV,  Non-Ideal  Detector  Faces,  addresses  both  of  these 
assumptions,  and  provides  several  simulations  to  examine  how  well  the  this  detector 
can  detect  light. 

Chapter  V,  Measurement  Theory,  examines  how  this  detector  can  be  used  in  the 
IGAR  system  to  measure  the  DHR  of  a  sample.  This  chapter  serves  as  the  bridge 
between  the  theory  of  the  operation  of  IGAR  and  implementing  this  system  in  a 
laboratory  setting.  Actually  implementing  this  system  is  covered  in  Chapter  VI, 
Equipment  Description.  Chapter  VII  two  of  the  most  important  pieces  of  software 
used  in  the  previous  chapter.  Finally,  Chapter  VIII,  examines  some  of  the  mea¬ 
surements  made  by  this  system,  and  compares  them  conceptually  to  the  simulations 
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performed  earlier. 


1.2  Axis  and  Units 

The  axis  shown  in  Figure  1  will  be  used  throughout  the  rest  of  the  document.  9 
refers  to  the  polar  angle  between  a  vector  and  the  z  axis.  0  is  the  azimuthal  angle, 
measured  from  the  x  axis  to  the  projection  of  the  vector  into  the  xy  plane.  For 
example,  a  vector  with  6  =  0°  would  be  parallel  to  the  normal  vector,  and  a  vector 
with  9  =  90°  would  be  perpendicular  to  the  normal. 
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Figure  1.  Axis  used  for  this  paper. 


The  SI  units  of  measure  are  used  throughout  the  entire  document.  Most  distances 
measured  are  relatively  short,  making  millimeters  the  preferred  unit  of  distance.  This 
makes  millimeters-squared  the  natural  unit  for  area.  Angles  will  be  defined  in  either 
radians  or  degrees.  Radians  are  the  preferred  unit  of  measure  in  Chapter  111,  due  to 
the  mathematics  involved  in  the  calculations.  Throughout  the  rest  of  the  document, 
degrees  are  used. 

All  power  measurements  are  reported  in  watts.  Wavelengths  will  be  measured  in 
microns.  DHR  measurements  are  unitless,  and  BRDF  measurements  have  units  of 
inverse  steradians. 
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II.  Discussion  of  DHR 


Directional  Hemispherical  Reflectance  (DHR)  is  a  measure  of  the  ratio  between 
the  light  incident  and  the  light  reflected  off  a  surface  for  a  single  angle  of  incidence 
[18].  Nicodemus  defined  this  relationship  as  the  flux  seen  looking  at  the  surface  over 
the  flux  incident: 

P(a;^)=d$ ?  (1) 

where  p  is  the  DHR  of  the  surface  in  question,  is  a  solid  angle,  the  subscripts  i 
and  r  refer  to  incident  and  reflected  radiation  respectively,  and  is  the  radiant  flux. 
This  is  depicted  in  Figure  2.  In  practice,  flr  typically  refers  to  the  entire  hemisphere 
of  reflectance,  making  DHR  a  function  of  only  angle  of  incidence.  The  majority  of 
samples  are  isotropic,  meaning  that  this  measure  can  be  defined  only  as  a  function  of 
angle  of  incidence  in  a  single  plane.  Samples  without  this  symmetry  must  be  defined 
for  both  this  angle  of  incidence  and  the  orientation  of  the  sample.  It  is  important  to 
note  that  several  other  variables  can  be  added  to  this  equation,  taking  into  account 
such  effects  as  sample  temperature,  input  polarization  or  wavelength.  These  effects 
will  be  ignored  throughout  the  paper,  and  DHR  will  only  be  talked  about  as  a  function 
of  angle  of  incidence  in  a  single  plane. 

DHR  can  be  expressed  in  terms  of  Bidirectional  Reflectance  Distribution  Function 
(BRDF)  [18], 


p{pjp  2n) 


J'u,i 


(2) 


where  fr  is  the  BRDF  of  the  material.  According  to  Nicodemus  [18],  if  we  assume  a 
uniform  and  isotropic  incident  beam,  and  take  the  wavelength  into  account,  we  can 
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Figure  2.  Picture  of  DHR  from  Nicodemus  paper.  Note  that,  in  the  rest  of  this  paper, 
DHR  refers  to  the  light  reflected  over  the  entire  hemisphere,  not  some  portion  uir. 

rewrite  the  equation  as 

P7t/2  /»27T 

p(6i,  4>i]  X)  =  /  /  fr(0r,  4>r]  9i4>i]  A)  cos6r  sindrd(j)rddr.  (3) 

Jo  Jo 

Note  that  while  this  equation  shows  a  clear  path  between  BRDF  and  DHR,  this 
calculation  is  non-trivial  given  real  data.  The  majority  of  measured  BRDF  data  is 
measured  for  only  a  single  value  of  <pr.  In  order  to  actually  calculate  the  integral,  the 
value  of  the  BRDF  for  all  other  values  of  6r  have  to  be  calculated.  Circular  symmetry 
around  the  incident  beam  is  a  good  assumption,  but  very  difficult  to  calculate  due 
to  the  geometry  involved.  There  are  very  few  practical  ways  to  go  between  measured 
BRDF  and  DHR  data.  Balling  did  make  this  conversion  in  his  research  [2],  but  only 
for  =  0.  The  angle  of  incidence  he  chose  made  the  geometry  of  the  problem  much 
easier. 

As  a  note,  DHR  and  Hemispherical  Directional  Reflectance  (HDR)  are  the  same 
information  measured  in  different  ways  [13].  In  a  DHR  measurement,  the  sample 
is  illuminated  from  a  single  angle,  and  the  signal  is  integrated  over  all  angles.  In 
HDR,  the  sample  is  illuminated  from  all  directions,  typically  using  a  blackbody  and 
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hemi-ellipsoid  mirror,  and  the  signal  is  measured  from  a  single  angle.  The  SOC-lOO, 
manufactured  by  Surface  Optics  Corp.,  San  Diego,  CA,  is  an  example  of  an  instrument 
which  operates  in  this  manner.  By  reciprocity,  both  of  these  measurements  yield 
the  same  information.  One  of  the  biggest  advantages  of  this  piece  of  equipment 
is  the  ability  to  measure  the  spectral  characteristics  of  the  sample.  Because  the 
sample  is  illuminated  by  a  blackbody,  the  instrument  can  measure  a  broad  range 
of  wavelengths  simultaneously.  However,  the  directional  emittance  of  the  blackbody 
limits  the  angular  capabilities  of  the  measurement,  necessitating  the  use  of  laser 
sources  in  a  grazing  angle  instrument. 
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III.  Detector  Theory 


The  purpose  of  this  section  is  to  prove  that  the  total  power  emitted  by  the  sample 
can  be  measured  by  this  detector.  In  this  experiment,  a  laser  is  used  to  illuminate  a 
sample,  which  is  located  at  the  bottom  focus  of  a  hemi-ellipsoidal  mirror.  The  light 
reflected  by  the  sample  is  scattered  onto  an  mirror,  and  reflected  to  the  second  focus 
where  a  detector  has  been  mounted.  This  detector  is  particularly  novel,  as  it  is  a  five 
sided  pyroelectric  detector.  The  goal  of  this  chapter  is  to  show  that  the  total  power 
emitted  by  the  sample  is  equal  to  the  sum  of  the  power  detected  by  each  side  of  the 
detector. 

Initially,  two  critical  assumptions  are  made.  The  first  assumption  is  that  the 
distance  from  the  sample  to  the  detector  doesn’t  change  over  the  size  of  the  of  the 
sample.  This  greatly  simplifies  the  problem  and  allows  a  closed  form  solution  to  the 
calculation.  In  general,  this  is  not  true,  and  will  be  addressed  in  Section  ??.  The 
second  assumption  is  that  the  detector  is  a  perfect  detector  and  can  perfectly  detect 
all  light  incident  upon  it.  This  is  also  not  true,  and  can  easily  be  addressed  after 
redefining  the  problem  to  address  the  previous  assumption.  By  making  these  two 
assumptions,  a  closed  form  solution  to  the  problem  can  be  calculated  and  compared 
to  the  ideal  case,  showing  that,  in  the  best  case  scenario,  this  detector  can  properly 
measure  the  power  reflected  by  the  sample. 

3.1  Sample  Power  at  the  Detector 

In  order  to  prove  that  the  total  power  reflected  by  the  sample  can  be  calculated 
by  summing  the  power  seen  by  each  side  of  the  detector,  the  power  seen  by  each 
side  of  the  detector  must  be  calculated.  Using  the  relationship  between  power  and 
radiance,  the  power  incident  on  each  side  of  the  detector  can  easily  be  calculated 
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given  any  particular  scatter  pattern.  Assuming  that  the  total  power  is  equal  to  the 
power  collected  by  each  face  of  the  detector,  the  sum  of  all  five  faces  should  be  equal 
to  the  total  power  emitted  by  the  sample. 

As  stated  before,  the  experiment  is  constructed  such  that  the  sample  and  detector 
are  placed  at  opposite  foci  of  an  hemi-ellipsoidal  mirror.  An  ellipsoidal  mirror  is 
constructed  in  such  a  way  that  any  ray  leaving  the  first  focus  will  be  reflected  to  the 
second  focus  [10].  This  means  that  the  detector  is  essentially  imaged  at  the  sample. 
Every  ray  leaving  the  sample  will  be  reflected  from  the  mirror  towards  the  detector. 
A  depiction  of  the  imaged  detector  can  be  seen  in  Figure  3. 

It  is  important  to  note  that  this  theoretical  description  of  the  system  assumes  that, 
from  the  point  of  view  of  these  radiometric  calculations,  the  image  of  the  detector  is  a 
perfect  image.  In  practice,  this  is  not  true.  The  hemi-ellipsoidal  mirror  deforms  each 
face  of  the  detector  into  a  curved  shape.  This  new  shape  is  very  difficult  to  analyze 
analytically,  necessitating  the  use  of  the  perfect  image  assumption.  The  small  size  of 
the  detector  relative  to  the  radii  of  the  hemi-ellipsoidal  mirror  make  this  deformation 
minor  and  this  assumption  holds  true.  Our  future  research  is  going  to  examine  the 
effects  of  accounting  for  the  deformation  of  the  image  of  the  detector  caused  by  the 
mirror.  With  the  deformed  image,  the  same  analysis  can  then  be  applied  to  the 
system. 

Because  of  the  setup  of  this  experiment,  the  system  can  be  approximated  as  an 
emitter  which  emits  light  in  the  same  way  that  the  sample  would  reflect  it.  The  mirror 
will  create  a  1:1  image  of  the  detector  in  object  space,  which  means  the  system  can 
be  modeled  as  if  the  detector  was  built  around  the  sample.  Note  that  this  approach 
ignores  complications  such  as  mirror  aberrations  and  imperfections.  To  account  for 
these  effects,  an  optics  modeling  package  such  as  Zernax  would  be  needed.  The  image 
of  the  detector  with  the  associated  dimensions  can  be  seen  in  Figure  4. 


Figure  3.  The  ellipsoidal  mirror  essentially  images  the  detector  into  object  space  around 
the  sample.  Assuming  the  sample  is  very  small,  the  detector  can  be  analyzed  as  if  it  is 
being  illuminated  by  a  point  source  located  at  the  foci  of  the  mirror. 


Figure  4.  This  is  a  diagram  of  the  image  of  the  detector.  The  face  perpendicular  to 
the  x  axis  is  the  Left  Face,  the  one  perpendicular  to  the  y  axis  is  the  Front  Face,  and 
the  one  perpendicular  to  the  2  axis  is  the  Top  Face. 


This  system  can  be  modeled  as  a  system  with  a  very  small  sample  located  at  the 
focal  point  of  the  mirror.  Specifically,  the  sample  must  be  small  enough  that  the 
distance  from  the  sample  to  the  detector  doesn’t  change  over  the  area  of  the  sample. 
Initially,  a  Lambertian  scatterer  will  be  analyzed.  This  scatterer  has  the  simplest 
definition,  simplifying  the  problem  and  enabling  a  closed  form  solution.  Later,  a 
simulation  will  be  used  to  analyze  different  scatter  patterns.  For  this  sample,  the 
energy  radiated  from  it  can  be  calculated  from  the  equation 

<92<L  =  L0  cos  6d  dAs  d VLd  (4) 
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where  <f>  is  the  total  power,  L0  is  some  constant  radiance,  As  is  the  area  of  the  emitter, 
and  fid,  is  the  solid  angle  in  which  power  is  emitted  [8].  To  calculate  the  total  power 
emitted,  we  need  to  integrate  over  both  the  area  of  the  emitter  and  the  solid  angle. 
Integrating  the  above  equation  over  the  area  of  the  detector  yields 

<9$  =  As  (  L0  cos  6  dQd.  (5) 

J  Q,s 

dnd  can  be  defined  in  terms  of  (j)d  and  9S: 


dttd  =  sin  9 d  d(j)d  ddd, 


(6) 


which  makes  (5) 

/*27T  f*7V  /  2 

® Optimal  =  AS  /  L0  cos  4>dsm (j)d  d,9d  d(j)d  =  n ASL0.  (7) 

Jo  Jo 

This  result  is  the  power  emitted  by  a  perfectly  Lambertian  scatterer. 

Now,  in  order  to  calculate  the  total  power  detected  by  the  detector,  the  power 
seen  by  each  face  must  be  calculated.  (5)  is  an  ideal  starting  point,  but  integrating 
over  angle  space  is  difficult.  Instead  a  rectangular  face  located  a  distance  away  from 
the  detector  can  be  used  instead  of  integrating  over  all  possible  angles,  as  seen  in 
Figure  5. 

Now,  dnd  will  be  redefined  as 


dnd  = 


dxd  dyd 


r 


2 


cos  a 


(8) 


where  r  is  the  distance  from  the  center  of  the  sample  to  the  center  of  the  differential 
area,  cos  a  is  the  angle  between  the  normal  of  the  differential  area  and  the  ray  going 
from  the  differential  area  to  the  sample.  Since  this  face  is  parallel  to  the  sample,  the 
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I/Max 


dy 

dx 


rW  z 


Figure  5.  The  top  face  is  essentially  imaged  upside  down  into  object  space.  An  ellip¬ 
soidal  mirror  makes  a  1:1  image  of  the  object,  meaning  that  the  image  of  the  detector 
has  the  same  dimensions  as  the  object. 


angle  a  is  the  same  as  the  angle  6d-  r  is  the  distance  from  onr  sample  (located  at 


(0,0,0)),  which  makes  the  differential  solid  angle 


_  dxd  dyd 

oild  —  — o  i — 2  cos  V d 

xd  +  Vd  +  4 


Since  this  face  faces  —z,  this  cosine  term  becomes  -.  Now,  instead  of  integrating  over 
angle  space,  we  can  integrate  over  ( Xd ,  Vd,  zd )  space. 


‘h'T'nn  — "  As 


cy  Max  r^Max 


'  2/Min  J  ^Min 


dxd  dyd 

L0  cos  0d^— — y- — o  cos  Vd- 

Xd  +  Vd  +  4 


We  choose  to  convert  cos  6d  to  cartesian  space  to  simplify  the  integration. 


a  " “ 

cos  0d  =  —  = 


'4  +  y2d  +  4 


which  makes  the  previous  equation 


!•  J/Max  rx  Max  J2 

$ToP  =  AS  /  L0  ===-  dx  dy. 

JyMin  dxMin  (  \J  X1  +  |/2  +  Z1 ) 2 
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For  simplicity,  a  detector  symmetric  about  the  axis  is  assume,  making  ^Min  =  — a^Max 
and  j/Min  =  —  2/Max-  z  is  a,  constant  over  the  surface,  and  is  equal  to  ZMax-  Completing 
the  integration  reveals 


®Top  —  A «L 


s-^0 


2^Max 


V^vr-v  + 


:  tan 


Max  '  ^Max 

2?/Max 


VVMax  +  Z2 


tan 


-l 


•^Max 


Max 


V  I/Max  +  " 


Max  . 


(13) 


Making  both  xMax  and  yMax  »  ^Max  gives 


‘f’Top  —  A,L 


S-OQ 


=  A,L 


s-OQ 


=  nAsL0 ,  (14) 


which  is  what  we  would  expect  for  a  Lambertian  source. 

Now,  for  the  sides  of  the  detector.  The  first  face  to  be  analyzed  will  be  the  front 
face  of  the  detector  (yd  is  held  constant).  This  is  depicted  in  Figure  6.  By  symmetry, 
the  back  is  the  same  as  the  front.  By  extension,  the  left  is  the  same  as  the  front  with 
^Max  and  I/Max  exchanged. 

(5)  is  again  used  as  a  starting  point.  The  same  transition  from  angular  space  to 
cartesian  space  is  used. 

dxd  dzd 

o\ld  — - - — cosa  (15) 

where  r  is  defined  the  same  as  above.  Because  the  face  points  in  the  —y  direction, 
cos  a  becomes  yd/r : 


Vd  _  _ Vd _ 

r  _  'M+W+4' 


(16) 
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Figure  6.  The  front  of  the  detector  is  imaged  in  the  same  fashion  as  the  top.  See  the 
caption  associated  with  Figure  5. 
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to  the  sum  of  power  seem  by  each  face  of  the  detector.  This  sum  of  powers  would  be 
equal  to  the  ideal  case,  which  is  7r AsL0. 

$2  =  <f>Top  +  d’Left  +  ^  Front  +  ^Right  +  ‘f’Back  (21) 

Due  to  the  symmetry  of  the  detector,  <l>Right  =  'I’Left  and  ^Back  =  ‘I’Front-  This  makes 
the  above  equation 

=  <f>Top  +  2<hLeft  +  2<hFront  (22) 

Inserting  the  equations  calculated  above  for  these  terms,  the  two  terms  of  <l>Top  cancel 
out  the  second  portion  of  both  2<hLeft  and  2<f>Front,  leaving 

=  ASL0  [2  tan-1  ( +  2tan"1  ( .  (23) 

\  ?/Max  /  \  ^Max  /  _ 

Defining  a  new  arbitrary  angle  (3  to  be  equal  to  tan~ 1  (xMax/i/Max)  •  H  follows  then 
that  7r / 2  —  /3  =  tan_1(|/Max/^Max),  making  the  above  equation 

=  ASL0  [2 f3  -  2/3  +  2vr/2]  =  nAsL0.  (24) 

In  (7),  this  is  shown  to  be  the  total  power  emitted  by  the  sample.  Therefore,  by 
summing  the  power  seen  on  each  side  of  the  detector,  the  detector  can  capture  the 
total  power  emitted  by  the  sample. 

3.2  Simulating  Large  Samples 

There  were  two  critical  assumptions  made  in  the  above  section.  The  first  is  that 
the  sample  is  significantly  smaller  than  the  detector.  This  allowed  the  formula  for 
calculating  the  distance  from  the  sample  to  the  detector  to  be  greatly  simplified, 
which  makes  the  analytic  solution  possible.  In  general,  this  is  not  the  case.  From 
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a  design  perspective,  the  detector  should  be  as  small  as  possible,  which  means  the 
radius  distance  will  vary  significantly  from  one  side  of  the  sample  to  the  other.  The 
other  assumption  made  was  that  the  detector  can  absorb  all  of  the  light  incident  upon 
it,  which  will  be  addressed  in  the  following  section. 

Both  of  these  assumptions  are  what  enabled  a  closed  form  solution  for  the  power 
seen  by  the  detector.  Now,  with  the  added  complexity  of  the  system,  a  numeric  solu¬ 
tion  is  the  preferred  solution  to  this  problem.  To  facilitate  this,  (4)  can  be  rewritten 
as 

A2<f>  =  L(9S,  4>s)  cos  9S AQdAAs  (25) 

Each  face  needs  to  be  treated  separately  from  the  other,  so  we  actually  need  to 
calculate  this  power  five  different  times.  We  will  use  the  subscript  j  to  indicate  which 
face  we  are  dealing  with.  A  list  of  the  index  j  and  its  associated  face  can  be  found  in 
Table  1.  This  makes  equation  (25) 

A2<f>j  =  L(9S,  4>s)  cos  9s/\£ldj^As.  (26) 


We  can  rewrite  A Qd,j  as 


Aftdj  = 


A  Adj 

r2(AAs,  A Ad) 


cos  9d- 


which  makes  (26) 


(27) 


A2<S  =  L($S)  cj>s)  cos  6S  cos  0d  1  A  1,/A/t  ,  (28) 

which  says  that  the  differential  power  is  a  function  of  both  the  position  of  the  differ¬ 
ential  area  of  the  sample  and  the  position  of  the  differential  area  of  the  detector. 

To  make  the  calculations  of  r(AAs,  A  Ad)  and  cos  6  easier,  we  will  now  define  the 
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Table  1.  List  of  the  index  j  and  it’s  associated  face. 


j  =1:  Front 
=  2:  Left 
=  3:  Back 
=  4:  Right 
=  5:  Top 


Figure  7.  The  vector  V  is  the  vector  starting  at  the  differential  area  of  the  sample  and 
ending  at  the  differential  area  of  the  detector. 

vector  V  as  the  vector  starting  at  the  differential  area  of  the  sample  and  ending  at 
the  differential  area  of  the  detector,  as  shown  in  Figure  7. 


V (A Ad,  A As)  =  V (xs,  ys ;  xd,  yd,  zd)  =  (. xd  -  xs)x  +  (yd  -  ys)y  +  (zd)z  (29) 

where  the  s  subscript  represents  the  position  on  the  sample,  and  the  d  subscript 
represents  the  position  on  the  detector.  Now,  the  distance  r  is  just  the  length  of  the 
vector  V.  (Note:  At  this  point,  the  (A Ad,  AAS)  portion  of  the  V  term  will  be  omitted 
for  brevity.) 


r(AAd,  AAS) 


\f  {xd  ~  xe)2  +  (yd  -  ye)2  +  zj,  (30) 
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The  term  cos  9S  can  now  be  written  as 


cos  9  s  = 


(31) 


cos  Ofij  now  becomes 

V  ■  n- 

cos  9dj  =  — =r-^-,  (32) 

"  ll^ll 

where  n  is  the  normal  vector  for  the  face  of  the  detector. 

In  order  to  calculate  the  power  incident  on  each  face  of  the  detector,  we  need  to 
apply  a  summation  to  (28). 


V-zV 


n-i 


As  A, 


d,j 


AAdJAAe 


(33) 


where  As  is  the  area  of  the  sample,  AAS  is  the  differential  area  of  the  sample,  Adj  is 
the  area  of  the  jth  detector  face,  and  A Adj  is  the  differential  of  that  area. 

As  is  the  range  of  x  values  from  —  £s,max  :  £S)inax  with  a  step  size  of  Axs,  and  the  y 
values  from  —  ys,max  ■  Us, max  with  a  step  size  of  A ys.  The  range  of  Adj  and  the  value 
of  A Ad  j  are  defined  in  Table  2.  This  formula  still  assumes  a  perfect  detector  which 
can  detect  all  of  the  light  incident  upon  it. 

As  shown  in  (21),  the  total  power  is  equal  to  the  sum  of  the  power  detected  by 
each  face  of  the  detector.  This  makes  the  total  power  of  the  detector 


$act  =  (34) 

j 

Performing  this  calculation  for  an  arbitrary  detector  and  sample  geometry  gives  re¬ 
sults  very  close  to  the  optimal  result  of  <f>det  =  7T AsL0 .  These  results  show  that  this 
detector  can  still  recover  all  of  the  power  emitted  by  the  sample,  which  is  what  is 
expected  given  the  results  from  the  previous  section. 
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Table  2.  List  of  the  index  j  and  the  information  associated  with  the  summation  for 
calculating  <!>, . 

FrOllt.  AdA  =  \jjd  =  ?/max  •  t/maxi  Zd  =  0  .  *max} 

Ai4d,i  =  AydAzd 
Left:  Ad)  2  =  {xd  = 

^max  •  ^max?  Zd  =  o  :  ^max} 

AAd)2  =  AxdAzd 
Back:  Ad )3  =  {yd  =  -y 

max  •  J/max)  Zd  =  0  :  ^max} 

AAd^  =  AydAzd 
Right:  AdA  =  { xd  = 

^max  •  ^max;  zd  =  o  •  ^max} 

AAdA  =  AxdAzd 

Top.  Ad,  5  {xd  2-max  •  ^maxj  Vd  I/max  •  I/max  J 

AAdA  =  AxdAyd 


%d  ^max 
h\  =  —  X 

Vd  y  max 

n2  =  -y 

%d  =  2- max 

fl'i  =  X 
yd  y  max 

h4  =  y 

cj  *max 

fir,  =  ~Z 


Now,  the  equation  (33)  addresses  the  first  assumption  made  in  the  previous  sec¬ 
tion.  This  equation  assumes  a  large  sample,  yet  shows  that  the  detector  can  still 
detect  all  of  the  light  scattered  by  the  sample.  This  new  framework  can  be  easily 
extended  to  address  the  second  assumption  that  was  made  earlier.  The  vector  V 
enables  us  to  implement  an  angle  dependent  absorption  function  on  each  face  of  the 
detector,  addressing  the  fact  that  no  real  detector  is  a  perfect  absorber. 
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IV.  Non-Ideal  Detector  Faces 


In  Chapter  III,  where  the  general  solution  for  the  power  seen  by  the  detector  was 
derived,  there  were  two  variables  which  were  undefined.  The  first  was  L(V),  which 
is  the  radiance  seen  by  the  differential  area  of  the  detector.  The  second  was  oi(V), 
which  is  the  absorptance  experienced  by  that  same  area.  The  radiance  can  be  defined 
very  well  using  a  BRDF  model,  and  the  absorptance  can  be  modeled  using  a  Fresnel 
absorptance.  After  defining  both  of  these  values,  we  can  calculate  the  power  seen  by 
the  detector,  and  analyze  it’s  behavior. 

4.1  Detector  Absorption 

The  second  assumption  that  was  previously  made  is  that  the  detector  can  absorb 
all  of  the  incident  light.  This  is  not  true.  While  dark  materials  absorb  most  energy, 
some  of  the  light  is  reflected.  This  implies  that  the  absorptivity  of  the  detector  is  a 
function  of  angle,  which  drastically  complicates  the  calculation  of  power.  Ideally,  this 
calculation  would  be  performed  with  a  BRDF  model  of  the  face  of  the  detector.  This 
paper  will  examine  the  case  of  a  flat  detector  with  some  index  of  refraction  n.  This 
air-detector  interaction  will  result  in  a  Fresnel  reflection,  from  which  the  absorption 
of  the  detector  can  be  calculated. 

The  Fresnel  reflection  of  the  detector  face  is  composed  of  two  terms:  one  for  the 
S  polarized  light  and  one  for  the  P  polarized  light  [4],  Both  of  these  terms  can  be 
calculated  as  a  function  of  the  angle  of  incidence  and  the  index  of  refraction  of  the 
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detector,  n.  (It  is  assumed  that  the  index  of  refraction  of  the  air  is  1.) 


Rs(n,6i) 


Rp(n,0i) 


cos  9i  — 

A- 

~2  sin2  9i 

nz  L 

cos  9i  +  n\ 

A- 

A  sin2  9~ 

nz  1 

1  —  -4  sin2  9j  —  n  cos  9, 


J  1—^2  shr2  9 i  +  n  cos  9i 


(35) 


(36) 


If  the  assumption  that  the  light  reflected  by  the  sample  is  randomly  polarized  is  made, 
the  reflectance  of  the  detector  face  will  be 


R(n,  9i) 


Rs(n,9j)  +  RP(n,9j ) 
2 


(37) 


Using  Kirchoff’s  Law  [14],  the  absorption  of  the  face  is  equal  to  one  minus  the  reflec¬ 
tion.  Therefore,  the  absorption  of  the  detector  face  can  be  expressed  as 


a(n,9i)  =  1  -  R(n,9i). 


(38) 


Note  that  the  above  equation  defines  a  in  terms  of  9%.  For  this  detector,  the 
angle  9^  is  equal  to  the  angle  between  the  differential  area  of  the  detector  and  the 
differential  are  of  the  sample,  which  is  9rjj .  Using  (38),  this  angle  can  be  calculated 
for  any  vector  V,  allowing  (38)  to  be  expressed  as 

a(V,  n)  —  1  —  R(V,  n).  (39) 

Now,  (33)  can  be  rewritten  to  account  for  the  absorptivity  of  the  detector. 

=  L(edJ^)a(V,  n) {V  '  "j)  A AjjAA..  (40) 
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Applying  the  same  summation  as  before, 


%  =  E  E  mj.mv,  n) (V'  •  %)  AAd, A  A.  (41) 

As  AdJ  1 1  ^  1 1 

If  a(V,n )  was  a  constant,  it  could  factored  out  of  the  summation,  and  could  be 
corrected  by  using  a  compensation  term.  However,  it  is  a  function  of  the  vector  V, 
which  means  that  information  is  lost,  and  cannot  correct.  Once  again,  the  total  power 
seen  by  the  detector  is  simply  the  sum  of  the  powers  seen  by  each  face  of  the  detector, 
or 

$det  =  £>;•  (42) 

3 

4.2  Samples 

Now,  the  only  thing  left  to  define  in  (41)  is  the  radiance  of  the  sample.  We  can 
easily  do  this  by  using  a  BRDF  model.  Radiance  is  related  to  BRDF  by  the  formula 

Lr{0,  fa  =  F(0i,  fa;  6,  flEiiOi,  fa)  cos  0i.  (43) 

where  Lr  is  the  reflected  radiance,  F  is  a  BRDF,  E,  is  the  incident  irradiance,  and  dj 
is  the  angle  of  the  incident  radiance  [3]. 

I  chose  to  model  four  samples  with  two  separate  BRDFs  for  this  simulation:  a 
perfectly  Lambertian  scatterer  and  three  samples  using  a  Cook-Torrance  model.  A 
Lambertian  scatterer  [12]  has  a  BRDF  of 

F(dl,fa;0,fa  =  -.  (44) 

7T 

While  this  model  is  heavily  simplistic,  it  gives  us  a  good  baseline  for  measuring  diffuse 
samples.  A  huge  number  of  diffuse  samples  can  be  approximated  as  Lambertian  or 
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Lambertian  like  scattering  patterns. 

The  Cook- Torrance  BRDF  model  is  a  significantly  more  complex  model  used  for 
simulating  specular  samples  [1]. 


K 


Ps 

7T 


DG 

(N-L)(N-V) 


Fresnel (F0,  V  ■  H) 


(45) 


where  N  is  the  normal  of  the  sample,  L  is  the  vector  of  the  incident  light,  V  is  the 
vector  from  the  sample  to  the  point  being  measured  (in  our  case,  from  the  differential 
area  of  the  detector  to  the  differential  area  of  the  sample),  H  is  the  half  angle  vector 
between  V  and  L ,  G  and  D  can  be  found  in  the  equations  below,  ps  is  one  of  the 
fitting  parameters,  and  Fresneli)  is  calculates  the  Fresnel  reflection  of  the  sample 
[15].  (For  a  more  complete  description  of  this  function  and  the  associated  variables, 
see  Cook’s  paper  “A  Reflectance  Model  for  Computer  Graphics”  [16].) 


G  =  min  <  1, 


D  = 


2 (N-H)(N-V)  2 (N-H)(N-L) 


(V-H) 


(V-H) 


m 2  cos4  5 


exp  — 


tan  5 


m 


(46) 

(47) 


The  D  term  contains  one  of  the  two  fitting  parameters  of  the  model,  the  variable  m. 
5  is  the  angle  between  the  half  angle  vector  and  the  normal.  “Experimental  Analysis 
of  BRDF  Models  -  Supplemental”  by  Ngan  et  al.  [1]  gives  examples  of  this  BRDF 
model  with  a  list  of  the  fitting  parameters  used.  Since  the  goal  is  to  test  the  detector 
and  not  the  sample  in  this  model,  several  samples  were  chosen  from  the  paper:  a 
diffuse  sample,  a  specular  sample,  and  one  sharing  both  characteristics. 

The  Fresnel  portion  of  the  equation  is  related  to  the  Fresnel  reflection  of  the 
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sample  [15],  and  is  calculated  using  the  following  formula: 


where 


Fresnel (F0,  V  ■  H) 


1  ( 9  ~  c)2  [1  ,  [c(9  +  c)  -  l]2' 

2  (g  +  cY[  +  [c(g-c)  +  l]2. 


(48) 


c  =  cos  9  =  V  -H 

(49) 

g2  —  n2  +  c2  —  1 

(50) 

1  +  \/Fq 

(51) 

n  =  - 

1  —  t/Fq 

Examining  these  equations  sheds  very  little  light  on  the  operation  of  the  shader,  and 
even  less  on  the  operation  of  the  detector,  so  1  will  not  go  into  further  detail  here. 
What  is  important  is  given  an  index  of  refraction  for  the  sample,  the  viewing  vector 
V  and  the  input  vector  L,  the  Fresnel  reflection  of  the  sample  can  be  calculated. 

For  this  simulation,  1  used  four  different  BRDF  patterns.  The  first  was  the  per¬ 
fectly  Lambertian  case.  A  Lambertian  scatterer  provides  a  good  baseline  for  overall 
detector  efficiency,  and  is  great  to  study  due  to  the  lack  of  any  angle  dependence. 

The  other  three  samples  were  all  taken  from  “Experimental  Analysis  of  BRDF 
Models  -  Supplemental”  [1].  All  three  use  the  Cook- Torrance  model  detailed  above, 
and  have  different  levels  of  diffusivity.  Their  fitting  values  can  be  found  on  Table  3. 
The  sample  “acrylic  blue”  is  the  most  specular  of  the  three,  and  has  a  tight  specular 
lobe  where  most  of  the  energy  is  concentrated,  “black  oxidized  steel”  is  a  much  more 
diffuse  sample,  and  is  a  good  example  of  a  sample  which  distributes  light  unevenly 
across  all  five  faces  of  the  detector,  “ch  ball  gold  metallic  2”  is  midway  between  the 
two,  and  provides  a  sample  with  attributes  of  each  of  the  previous  two  simultaneously. 
The  BRDFs  of  all  four  scatter  patterns  can  be  seen  in  Figure  ??. 
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Table  3.  List  of  sample  scatter  pattern  fitting  values.  All  values  were  taken  from  Ngan’s 
paper  ”  Experimental  Analysis  of  BRDF  Models  -  Supplemental”  [1]. 

Lambertian:  F(6,  </>)  =  ! 


Acrylic  Blue 
CH  Ball  Gold  Metallic  2 
Black  Oxidized  Steel 


F0  =  0.035 
F0  =  0.629 
F0  =  0.035 


m  =  0.037 
m  =  0.144 
m  =  0.19 


(c)  Black  Oxidized  Steel  Scatterer. 


(d)  Gold  Metallic  Scatterer. 


Figure  8.  Sample  scatter  patterns.  All  figures  have  units  sr  1.  All  patterns  were 
evaluated  with  9  =  45°  and  <j>  =  45°. 
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4.3  Interpretting  Results 


Ultimately,  the  goal  of  this  section  is  to  compare  the  power  measured  by  the 
actual  detector  to  that  of  the  ideal  detector.  The  ideal  detector  power  is  calculated 
as  a  function  of  angle  using  the  scatter  patterns  outlined  above  and  doesn’t  include 
the  absorption  term  of  the  detector.  This  is  the  actual  power  reflected  by  the  sample, 
and  the  baseline  of  the  comparison.  The  actual  power  is  the  power  detected  by  the 
detector  taking  into  account  the  absorption  and  geometry  of  the  detector.  Ideally, 
these  two  values  should  be  identical. 

Because  of  the  absorption  term,  the  detector  must  be  calibrated  to  calculate  the 
DHR.  This  calibration  will  be  further  discussed  in  the  section  relating  to  Measurement 
Theory,  V.  For  now,  we  are  going  to  implement  a  single  scalar  calibration  term  for 
calculating  the  DHR  from  the  power  calculated  by  42  and  34.  These  calibration  terms 
are  as  follows: 


Cdet  —  <hdet(0,  0)  1 


(52) 


Cac,  =  *«t(0,0)  \ 


(53) 


where  the  gold  metallic  sample  is  used  as  the  scatterer.  Basically,  this  is  simply 
setting  the  measured  DHR  of  both  the  simulated  detector  DHR  and  the  actual  DHR 
equal  to  1  at  9  —  0°.  This  is  equivalent  to  what  is  done  later  on  the  Measurement 
Theory  section. 

4.4  Simulation  Parameters 

In  order  to  simulate  the  detector,  the  various  values  found  in  Table  2  have  to  be 
defined.  The  values  used  in  this  simulation  can  be  found  in  Table  4.  The  dimensions 
of  the  detector  were  taken  from  the  actual  experiment.  The  IGAR  detector  is  10mm 
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Table  4.  Relevant  values  for  the  simulation. 


■t'max 
I/max 
*max 
■t’.s.rna.x 
Vs,  max 


10 

Axd  = 

0.1 

20 

Ayd  = 

0.1 

10 

A  zd  = 

0.1 

0.5 

Axs  = 

0.1 

0.5 

Ays  = 

0.1 

wide  by  20mm  long  by  10mm  tall.  The  range  for  the  sample  was  chosen  such  that  the 
specular  reflection  of  the  sample  was  small  enough  to  display  the  angular  dependence 
of  the  detector.  The  A  values  were  chosen  to  maintain  accuracy  while  minimizing 
calculation  times.  Doubling  the  A  values  drastically  increased  the  time  necessary  for 
computations  while  having  negligible  effect  on  the  overall  calculations. 

Note  that  all  of  the  code  used  for  these  calculations  can  be  found  below  in  Ap¬ 
pendix  B. 


4.5  Results 


The  simplest  scatter  pattern  to  analyze  is  the  Lambertian  case.  For  this,  the 
actual  DHR  of  the  sample  (utilizing  the  calibration  terms)  is  1.0602.  This  error 
shows  that  the  gold  scatterer  used  for  calibration  does  not  have  a  DHR  of  1.  This 
kind  of  error  can  easily  be  removed  in  a  real  system  by  accurately  measuring  the  DHR 
of  the  calibration  mirror.  For  the  detected  DHR  (the  term  including  the  absorption 
of  the  detector),  the  measured  DHR  has  a  value  of  1.0228.  Defining  the  error  as 


$act  ~  $det 
^act  +  ^det 


(54) 


the  error  between  these  two  measurements  is  a  little  less  than  2%. 

This  result  shows  two  important  characteristics  of  this  detector  which  will  be  rein¬ 
forced  by  the  other  samples.  One,  the  DHR  measurement  is  not  perfect.  Even  for  the 
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simplest  scatterer,  there  is  error  associated  with  this  measurement  before  accounting 
for  other  sources  of  error.  This  error  is  caused  by  the  fact  that  the  reflectance,  and 
hence  absorptance,  of  the  detector  is  a  function  of  angle.  The  light  scattered  to  the 
edges  of  the  detector  are  not  detected  as  well  as  the  light  scattered  to  the  center  of 
the  face. 

The  second  important  characteristic  is  that  this  result  is  still  very  close.  Even 
with  poorly  calibrated  data,  the  results  are  still  within  2%  of  each  other.  Later,  it 
is  shown  that  this  error  is  significantly  less  than  the  error  caused  by  the  uncertainty 
of  the  measurement.  For  example,  the  error  in  measuring  the  voltage  of  the  detector 
accounts  for  uncertainties  as  great  as  4%  in  the  measurement  of  DHR. 

The  next  scatterer  to  examine  is  the  Black  Oxidized  Steel  scatterer.  This  is  a 
very  diffuse  sample,  so  results  similar  to  the  Lambertian  case  would  be  expected. 
The  actual  DHR  and  detected  DHR  can  both  be  found  in  Figure  9.  Here,  we  find 
that  both  of  these  values  are  very  close  to  each  other  throughout  the  entire  range  of 
6  from  1°  to  90°.  There  is  still  a  difference  between  the  two  measurements,  but  it  is 
still  very  close.  The  measurement  has  an  average  error  of  2.5%. 

The  Acrylic  Blue  scatter  has  a  much  tighter  specular  lobe  than  the  Black  Oxidized 
Steel  scatterer,  and  because  of  it  has  a  different  error,  as  seen  in  Figure  10.  The 
Black  Oxidized  Steel  scatterer  error  was  almost  completely  positive,  indicating  that 
the  actual  DHR  was  greater  than  the  detected  DHR.  The  opposite  is  true  for  the 
Acrylic  Blue  scatterer.  Throughout  the  entire  range  of  angles,  the  detected  DHR  was 
greater  than  the  actual  DHR.  This  scatterer  had  a  smaller  error  than  the  previous 
case,  with  an  average  error  of  1.7%. 

The  final  scatterer  analyzed  was  the  Metallic  Gold  scatterer.  This  scatterer  had 
the  most  unique  DHR  of  the  group  of  samples,  as  seen  in  Figure  11.  Even  with 
the  irregularity  of  the  DHR,  the  detector  was  still  able  to  make  a  very  accurate 
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Figure  9.  Actual  and  Detected  DHR  for  the  Black  Oxidized  Steel  scatterer.  The  solid 
line  is  the  detected  DHR,  while  the  dashed  line  is  the  actual  DHR  of  the  sample. 
ADHR  =  Actual  DHR  -  Detected  DHR.  The  detected  DHR  varies  from  less  than  the 
actual  DHR  to  greater  than  the  actual  DHR,  but  is  still  very  close. 
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Figure  10.  Actual  and  Detected  DHR  for  the  Acrylic  Blue  scatterer.  The  solid  line 
is  the  detected  DHR,  while  the  dashed  line  is  the  actual  DHR  of  the  sample.  ADHR 
=  Actual  DHR  -  Detected  DHR.  The  detected  DHR  is  always  slightly  less  than  the 
actual  DHR,  but  is  very  close  over  the  entire  range  of  angles. 
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Figure  11.  Actual  and  Detected  DHR  for  the  Gold  Metallic  scatterer.  The  solid  line 
is  the  detected  DHR,  while  the  dashed  line  is  the  actual  DHR  of  the  sample.  ADHR 
=  Actual  DHR  -  Detected  DHR.  Even  with  the  complicated  nature  of  this  DHR,  the 
Detected  and  Actual  DHR  are  still  very  close  to  each  other. 


measurement.  The  average  error  of  this  sample  was  2.5%.  The  most  interesting 
feature  of  this  sample  is  how  much  the  error  varied  over  the  range  of  angles  sampled. 
The  error  at  40°  is  noticeably  more  than  the  error  at  10°. 

4.6  Summary 

These  results  demonstrate  many  important  characteristics  of  the  performance  of 
the  detector.  Firstly,  the  detector  cannot  make  perfect  measurements.  Even  with  the 
sample  which  was  calibrated  to  itself,  there  was  a  nontrivial  error  associated  with  the 
measurement,  before  accounting  for  any  uncertainties.  However,  these  errors  are  low 
enough  that  they  can  be  ignored  in  all  situations. 

The  second  result  is  that  the  error  is  a  function  of  angle,  and  hence  cannot  be 
calibrated  or  corrected  out  of  the  system.  The  detector  measured  the  Blue  Acrylic 
scatterer  as  having  a  DHR  greater  than  the  actual  DHR,  while  the  Black  Oxidized 
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Metal  scatterer  was  measured  with  a  DHR  less  than  the  actual  DHR.  The  Gold 


Metallic  scatterer  has  a  DHR  significantly  different  from  the  other  two  analyzed 
scatterers,  and  hence  the  most  unique  error.  This  error  varied  significantly  as  a 
function  of  angle.  Due  to  the  variety  of  these  errors,  they  cannot  all  be  calibrated 
out  of  the  system. 

However,  in  all  cases,  the  measured  error  was  low.  The  highest  measured  error 
over  all  four  samples  was  2.55%.  In  practice,  this  error  is  negligible  compared  to  the 
uncertainty  in  the  measurements  taken  by  the  system.  Even  without  being  able  to 
absorb  all  of  the  light  emitted  by  the  sample,  the  theoretical  detector  can  still  make 
accurate  measurements  of  the  power  reflected  by  the  sample,  and  hence  the  DHR. 
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V.  Measurement  Theory 


Given  the  detector  described  in  ,  how  can  DHR  be  measured  in  a  laboratory  set¬ 
ting?  As  shown  in  equation  (1),  two  measurements  are  required:  the  power  reflected 
by  the  sample  and  the  power  incident  on  the  sample.  From  Section  3.1,  the  total 
power  reflected  by  the  sample  can  be  calculated  by  summing  the  power  seen  by  the 
detector,  but  this  glosses  over  an  important  aspect  of  detector  operation.  The  detec¬ 
tor  does  not  report  power,  it  reports  a  voltage.  We  need  to  carefully  calibrate  the 
detector  in  order  to  recover  the  power.  This  calibration  method  is  covered  in  the 
following  sections. 

In  the  sections  that  follow,  it  is  assumed  that  the  system  has  access  to  the  following 
measurements:  detector  voltage  as  a  function  of  angle  and  of  face  with  uncertainties 
( Vd,j(& )),  power  seen  at  the  wavemeter  with  uncertainties  (Pw),  and  the  reflectivity 
of  the  some  calibration  standard  (for  this  experiment,  calibration  mirrors  were  used) 
with  uncertainties  ( pc{6 )).  No  other  physical  parameters  of  the  system  need  to  be 
known  in  order  to  make  DHR  measurements  of  the  sample. 

5.1  Basic  DHR  Measurements 

The  method  that  IGAR  uses  to  calculate  DHR  is  derived  from  equation  (1). 
Knowing  the  power  incident  on  the  sample  and  the  power  reflected  by  the  sample, 
the  DHR  can  be  calculated.  This  can  be  stated  as 

Mg) 

where  p{9)  is  the  DHR  as  a  function  of  angle  of  incidence,  is  the  power  reflected 
off  the  sample,  and  is  the  power  incident  on  the  sample.  Neither  $r  nor  <f>j  can  be 
measured  directly.  Instead,  we  measure  the  voltage  at  the  detector  and  the  power  at 
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the  wavemeter.  We  can  describe  the  relationships  between  these  values  and  and 
as 


M0)  =  cdvd{0)  (56) 

=  CWPW  (57) 

where  Cd  and  Cw  are  linear  loss  and  conversion  terms,  Vd  is  the  voltage  measured  at 
the  detector,  and  Pw  is  the  power  seen  by  the  wavemeter.  (Note  that  Vd  is  actually  a 
five  element  vector,  one  voltage  per  face.  This  will  be  addressed  in  Section  5.5.)  Cd 
and  Cw  are  coefficients  that  encompass  many  linear  terms.  For  example,  Cw  includes 
the  losses  due  to  the  beamsplitter  and  power  lost  due  to  the  input  aperture  of  the 
wavemeter.  While  the  values  of  these  terms  can  be  calculated  explicitly,  it  is  easier 
to  combine  these  terms  and  rewrite  them  in  terms  of  measurable  data. 

5.2  Calculating  Cd  and  Cw 

Values  for  both  Cd  and  Cw  can  be  calculated  explicitly.  Both  of  these  constants 
are  related  to  the  optics,  efficiencies,  and  the  responsivity  of  the  system.  Cw  relates 
the  power  seen  by  the  wavemeter  to  the  power  output  by  the  laser. 

Cw  =  N\^Mp\Mp2B  (58) 

where  W.o  is  the  transmittance  of  the  ND  filter,  Mp\  and  Mp 2  are  the  reflectivities  of 
the  periscope  mirrors,  and  B  is  the  reflectance  of  the  beamsplitter.  Unless  the  laser 
power  is  high  enough  to  cause  damage  to  the  optics,  all  of  these  values  are  constants. 
This  means  that  (57)  is  always  linear. 

Calculating  the  value  of  Cd  is  slightly  more  difficult,  but  follows  the  same  pattern. 
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This  value  can  be  expressed  as 


Cd  —  RdMh(l)NrTx/2Tc(l  —  B)  (59) 

where  Rd  is  the  responsivity  of  the  detector,  Mh  is  the  reflectance  of  the  hemi- 
cllipsoidal  mirror,  Nr  is  the  transmittance  of  the  variable  ND  filter,  Tw 2  is  the  trans¬ 
mittance  of  the  half-wave  plate,  Tc  is  the  effective  transmittance  of  the  chopper,  and 
B  is  once  again  the  reflectivity  of  the  beamsplitter. 

The  only  portion  of  this  term  that  is  not  explicitly  linear  is  the  detector  respon¬ 
sivity.  The  assumption  that  this  value  is  linear  is  crucial  to  the  operation  of  1GAR 
as  it  currently  functions.  If  the  term  is  linear,  then  the  techniques  detailed  in  Section 
6.4  can  be  used  to  calculate  a  calibration  constant,  and  we  can  avoid  calculating  these 
values  entirely.  However,  if  the  term  is  not  linear,  there  is  a  non-linear  term  that  needs 
to  be  accounted  for.  Some  complete  description  of  the  voltage-to-power  relationship 
would  be  required  in  order  to  calculate  the  DHR  of  the  sample.  By  operating  the 
detector  with  low  incident  power,  the  assumption  can  be  made  that  the  responsivity 
of  the  detector  is  linear. 

5.3  Atmospheric  Effects 

One  subject  that  was  not  addressed  in  the  previous  section  was  the  effects  of  the 
atmosphere  on  the  system.  In  the  portion  of  the  spectrum  covered  by  this  experiment, 
the  atmosphere  has  a  high  amount  of  attenuation,  and  can  have  strong  local  effects 
[9].  For  example,  a  drop  in  pressure  overnight  can  lower  the  attenuation  of  the 
atmosphere.  Some  of  these  effects  can  be  calibrated  out  of  the  system,  while  others 
cannot. 

The  beam  train  can  experience  atmospheric  effect  in  three  different  places:  be- 
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fore  the  beam  splitter,  between  the  beam  splitter  and  the  power  measurement,  and 
between  the  beamsplitter  and  the  pyroelectric  detector.  If  the  attenuation  is  con¬ 
stant  over  time,  the  drop  in  power  can  be  calibrated  out  of  the  system.  The  extra 
attenuation  would  show  up  as  an  extra  term  in  either  Cw  or  Cd ,  or  would  simply 
look  like  a  drop  in  power  output  from  the  laser.  Recalibrating  the  system  with  the 
atmospheric  effects  would  remedy  the  problem  and  remove  this  source  of  error  from 
the  measurement. 

Attenuation  as  a  function  of  time  cannot  be  so  easily  be  corrected.  If  the  at¬ 
tenuation  occurs  before  the  beam  splitter,  the  system  can  still  make  an  accurate 
measurement.  The  system  will  observe  a  drop  in  power,  but  the  ratio  between  the 
wavemeter  and  the  pyroelectric  detector  will  remain  unchanged.  If  the  attenuation 
occurs  after  the  beam  splitter,  one  of  the  detectors  will  see  less  power  than  the  other, 
and  the  ratio  between  the  two  will  differ. 

In  practice,  these  problems  have  not  been  observed.  Data  collected  from  the  same 
sample  with  different  calibrations  remained  accurate  despite  possible  changes  in  the 
atmosphere.  Overall,  the  lack  of  change  in  our  measurement  proves  the  robustness  of 
our  calibration  methodology. 

Other  DHR  systems,  such  as  the  SOC-lOO,  circumvent  this  problem  by  forcing 
their  reference  beam  to  follow  nearly  the  same  path  as  their  beam  interrogating  the 
sample.  Because  both  beams  have  very  similar  paths,  they  are  both  attenuated  by 
the  same  amount.  The  sample  beam  is  normalized  by  the  reference  beam,  and  any 
attenuation  due  to  atmospheric  effects  is  canceled  out. 
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5.4  Calibration 


Using  equation  (56),  the  value  of  the  calibration  coefficient  can  be  expressed  as 


=  c(9c) 

'  Vdc(6c) ' 


(60) 


where  the  c  subscript  implies  a  value  measured  during  calibration.  Using  equation 
(55),  this  equation  becomes 

s-t  Pci@c)^ic 

d=  Vdc(9c)  ‘ 

equation  (57)  allows  this  term  to  be  rewrite  in  terms  of  measurable  data  and  a  single 
unknown: 


ri  Pc(0c)Pwc 

Ct=Cw^m-  (62) 

This  implies  that  the  value  of  Cd  can  be  calculated  in  terms  of  Cw  if  the  calibration 
DHR,  the  power  at  the  wavemeter  during  calibration,  and  the  voltage  at  the  detector 
are  known  during  the  calibration. 

Further  simplifications  will  allow  the  Cw  term  to  be  removed  from  the  equation 
entirely.  Substituting  (62)  into  (56), 


■tr(6)=Cwpc(Oc)^-Vd(0). 


Substituting  both  this  and  (57)  into  (55), 


m  =  Crcipc>  VPZM») 
*  dc\P c) 


which  simplifies  to 

MW* vm 

P[  ’  Kfc(ec)  pw 


(63) 


(64) 


(65) 
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Combining  the  calibration  terms  into  a  single  variable,  the  DHR  can  be  calculated  as 


Vd(9) 


where 

Cc  =  t.  (67) 

Vdc\vc ) 

What  this  means  is  that  IGAR  can  be  calibrated  only  using  data  measured  by 
the  system.  No  information  regarding  any  of  the  optics,  lasers,  or  detectors  need  to 
be  known  in  order  to  make  DHR  measurements.  This  affords  much  flexibility  for  the 
system.  The  laser  or  the  optics  of  the  system  can  be  changed  without  making  any 
changes  to  the  operation  of  the  system. 


5.5  Making  Measurements 

Equation  (66)  ignores  one  significant  aspect  of  IGAR.  There  is  not  a  single  detec¬ 
tor;  there  are  five.  However,  as  shown  in  (21),  the  total  power  reflected  by  the  sample 
is  equal  to  the  sum  of  the  power  detected  by  each  face  of  the  detector.  Unfortunately, 
these  powers  are  not  known.  The  only  data  recorded  by  the  system  is  the  voltage  of 
each  face  of  the  detector. 

Starting  again  at  equation  (56),  a  single  DHR  value  can  be  calculated  in  the 
following  manner.  In  this  equation,  <hr  is  the  power  seen  by  the  entire  detector.  The 
power  seen  by  the  entire  detector  is  simply  the  sum  of  all  the  power  seen  by  all  the 
faces.  This  can  be  written  as 

5  5 

$r(0)  =  Y,  ‘M0)  =  5Z  CdJyd,M)-  (68) 

3= 1  i=i 

In  this  equation,  $rj  is  the  power  seen  by  face  j,  Cdj  is  the  calibration  coefficient  of 
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face  j,  and  Vaj  is  the  voltage  seen  on  face  j. 

In  order  to  make  this  measurement,  we  need  to  calculate  Cdj.  This  can  be  done 
exactly  the  same  way  it  was  done  before.  Equation  (62)  can  be  rewritten  as  a  function 


of  face: 


Cdj  = 


cw 


P  c\@  C,j)  Pwc,j 
Vdc,j  (@c,j  ) 


(69) 


It  is  important  to  note  that  Cw  is  not  a  function  of  face,  as  it  is  derived  from  the 
optics  relating  to  the  wavemeter. 

Combining  (68),  (69),  and  (57)  into  (55), 


P(8) 


r  p 

w 


5 

£ 

j= i 


c. 


Pc.j  (9cj)Pt 


WC,J 


Vdc,j  \@c,j 


vd,M 


(70) 


The  Cw  terms  cancel,  and  all  the  calibration  terms  can  be  again  combined  into  a 
single  constant  as  a  function  of  face.  This  makes  the  equation  for  calculating  DHR 
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where 


Ccj  —  Pc,jific,j) 


p. 


WC,J 


Vdc,j(@c,j) 


(72) 


It  is  important  to  note  that  writing  the  calibration  constant  in  this  way  gives  complete 
control  over  how  the  machine  is  calibrated.  The  DHR  of  the  calibration  sample  and 
the  angle  of  incidence  can  both  be  chosen  as  a  function  of  face.  This  means  samples 
and  angles  can  be  chosen  such  that  the  voltage  on  the  face  are  maximized  during 
calibration,  giving  more  accurate  calibrations. 
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5.6  Uncertainty 


In  any  measurement,  there  is  some  uncertainty.  In  this  particular  measurement, 
there  are  five  different  sources  of  error:  voltage  measured  during  the  measurement, 
power  measured  during  the  measurement,  voltage  measured  during  calibration,  power 
measured  during  calibration,  and  the  DHR  of  the  calibration  mirror.  In  order  to  better 
understand  these  uncertainties  and  their  effect  on  the  measurement  of  the  DHR,  begin 
by  assuming  that  the  system  is  perfectly  calibrated.  Beginning  with  equation  (71), 
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which  expands  to 
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which  simplifies  to 
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Therefore,  the  error  in  DHR  for  the  system,  assuming  a  perfect  calibration,  is 
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A  more  complete  derivation  of  this  equation  can  be  found  in  Appendix  A. 

Calculating  the  error  for  the  calibration  coefficient  is  slightly  more  difficult.  In 
(72),  we  have  three  variables  with  uncertainties.  Acknowledging  uncertainties,  this 
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equation  becomes 


=  W  +  (77) 

The  j  designator  on  Cc,  pc,  and  Vdc  have  been  omitted  for  clarity.  Expanding  the 
division  term,  this  equation  becomes 
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Next,  we  need  to  multiply  out  the  p  +  A p  term: 
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Simplifying, 
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or 

ACc  =  77  7a  \  Pc  A  Vdc(9c)  +  pc(9c)APwc  +  P^cApc(6*c)] .  (81) 
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Now,  the  effects  of  uncertainty  in  calibration  on  the  measurement  can  be  exam¬ 
ined.  Going  back  to  equation  (76), 
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From  (81),  A Ccj  is  seen  to  be  a  term  composed  only  of  other  A  terms.  This  allows 
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the  equation  to  be  simplifcd  to 


=  m^  +  ±Cc,^. 

*  w  ,  -Lw 

3= 1 

Interestingly,  the  uncertainty  in  the  calibration  has  no  effect  on  the  uncertainty  of 
the  measurement  for  small  uncertainties.  (For  more  significant  uncertainties,  the 
ACcjAVdj  cross  term  cannot  be  ignored.  In  practice,  this  case  does  not  occur.) 

All  of  the  above  equations  are  only  true  if  the  uncertainties  are  much  smaller  than 
their  measurements.  For  Pw ,  Pwc,  pc,  and  Vdc ,  this  is  always  true.  The  laser  is  always 
strong  enough  to  give  the  wavemeter  a  good  measurement,  the  calibration  sample  is 
a  mirror  with  around  97%  reflectivity,  and  the  detector’s  voltage  during  calibration 
is  typically  fairly  high.  However,  this  assumption  is  not  always  true  to  the  voltage  on 
the  detector  face  during  sample  measurements.  Samples  with  low  DHR’s  have  low 
voltages  with  significant  uncertainties. 

5.7  IGAR  Algorithm 

Due  to  the  uncertainties  in  the  measurement,  simply  summing  the  powers  of  all  of 
the  faces  together  increases  the  uncertainty  of  the  DHR  measurement.  For  the  case 
of  a  diffuse  sample,  all  faces  give  the  system  information  about  the  sample,  as  all  five 
faces  can  see  light  reflected  by  the  sample.  This  is  not  the  case  for  all  samples.  For 
example,  a  specular  sample  might  illuminate  only  a  single  face.  In  this  case,  one  of 
the  faces  is  contributing  information  about  the  DHR  of  the  sample,  and  the  other 
four  are  contributing  noise. 

To  combat  this,  NGC  implemented  an  algorithm  to  minimize  this  uncertainty  [6]. 
This  algorithm  is  described  later  in  Section  7.2.  In  a  mathematical  sense,  it  functions 
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as  a  sifting  feature  inside  the  summation  of  equation  (71). 
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The  addition  of  this  function  inside  the  summation  makes  equation  (83)  into 
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If  the  face  in  question  has  no  light  incident  on  it,  the  voltage  will  be  close  to 
zero.  This  means  that  in  equation  (84),  it  will  be  contributing  almost  nothing  to 
the  DHR.  Therefore,  the  inclusion  of  the  Ij  term  will  have  no  effect  on  the  measured 
DHR.  However,  it  makes  a  significant  difference  in  equation  (85).  Even  with  no 
power  incident  on  a  face,  the  system  can  still  measure  uncertainty  for  that  face’s 
voltage.  The  I3  term  effectively  nullifies  that  portion  of  the  uncertainty,  meaning 
that  the  overall  uncertainty  of  the  measurement  is  decreased.  This  leads  to  overall 
more  accurate  DHR  measurements  by  intelligently  selecting  which  faces  are  measured 
for  any  particular  sample. 


5.8  Summary 

By  carefully  defining  how  IGAR  makes  measurements,  calibration  can  be  per¬ 
formed  without  any  previous  knowledge  of  the  optics  of  the  system.  This  allows 
maximum  flexibility  in  making  measurements.  The  laser  source  and  optics  can  easily 
be  adjusted  for  making  any  measurement  without  sacrificing  measurement  accuracy. 
This  measurement  theory  was  implemented  at  AFIT  and  used  to  measure  the  DHR 
of  a  number  of  samples. 
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VI.  Equipment  Description 


So  far,  only  the  theory  of  the  IGAR  system  has  been  discussed.  Actually  imple¬ 
menting  these  theories  in  equipment  is  a  challenging  task.  The  next  few  sections  will 
examine  some  of  the  pieces  of  equipment  used  to  realize  the  above  equations,  and 
address  the  problems  and  solutions  related  to  implementing  the  theories  presented  in 
this  paper. 

6.1  Theory  of  Operation 

The  goal  of  IGAR  is  simple.  We  want  to  calculate  the  DHR  of  a  sample  using 
equation  (84).  Looking  back  at  that  equation,  there  are  only  two  pieces  of  information 
which  are  required:  the  power  incident  on  the  sample  and  the  power  reflected  off 
the  sample.  This  can  be  accomplished  using  two  detectors.  A  Bristol  Wavemeter 
measures  laser  power  using  the  detector  used  in  measuring  the  wavelength.  A  novel 
five  sided  pyroelectric  detector  is  used  in  conjunction  with  a  hemiellipsoidal  mirror 
to  calculate  the  power  of  the  light  reflected  off  the  sample.  Additional  information  is 
collected  to  calculate  the  uncertainty  of  the  measurement,  so  that  the  equipment  can 
ultimately  say  this  is  the  DHR  of  the  sample  at  this  wavelength  and  angle  with  this 
amount  of  uncertainty. 

6.2  Hardware 

IGAR  can  be  decomposed  into  three  main  pieces:  the  tunable  laser,  the  hemi- 
cllipsoidal  mirror,  and  the  detector,  arranged  in  the  manner  shown  in  Figure  12.  A 
photo  of  IGAR  can  be  found  in  Figure  13.  Light  is  emitted  by  the  laser,  and  directly 
illuminates  the  sample  mounted  at  the  lower  focal  point  of  the  hemi-cllipsoidal  mirror, 
as  depicted  in  Figure  12.  A  portion  of  the  beam  is  diverted  away  for  measurements 
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Figure  12.  Functional  diagram  of  the  IGAR  system. 

performed  by  the  Bristol  741  Wavemeter,  manufactured  by  Bristol  Instruments  of 
Victor,  NY.  Both  the  mirror  and  the  sample  are  mounted  to  a  rotation  stage,  allowing 
control  of  the  angle  of  incidence  of  the  laser  light.  The  light  is  reflected  off  the  sample 
and  onto  the  mirror.  The  mirror  then  focuses  the  light  to  the  second  focal  point, 
where  the  pyroelectic  detector  is  mounted.  This  signal  is  measured  by  the  software, 
which  then  calculates  the  DHR  value. 


Figure  13.  Photograph  of  the  IGAR  system. 


Northrop  Grumman  chose  a  LASY-4G  tunable  CO2  laser,  manufactured  by  Access 
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Laser  of  Everett,  WA,  as  the  illumination  source  for  1GAR  [5].  This  laser  has  ten 
useable  laser  lines  ranging  from  9.271pm  to  10.718pm,  and  was  an  excellent  choice 
for  several  reasons.  First,  the  laser’s  tunability  allows  the  equipment  to  cover  a 
good  portion  of  the  LYV1R  atmospheric  window,  a  part  of  the  spectrum  of  particular 
interest  to  the  Air  Force.  Second,  the  spacing  of  the  laser  lines  allows  for  a  good 
spectral  coverage,  while  the  spacing  between  them  ensures  that  the  equipment  will 
not  accidentally  tune  to  the  wrong  wavelength.  The  low  power  of  the  laser  allows 
enough  light  to  illuminate  the  sample  while  never  damaging  the  material.  (The  spec 
sheet  lists  the  maximum  power  as  400mW,  while  in-house  measurement  places  this 
value  closer  to  40mW.  Individual  lines  can  be  even  weaker.)  Northrop  Grumman  did 
extensive  testing  regarding  the  stability  and  structure  of  each  wavelength’s  modes, 
ultimately  deciding  on  ten  wavelengths  with  the  best  structure.  The  laser  can  easily 
be  run  continuously,  and  the  software  allows  for  laser  tuning  automatically.  All  of 
this  means  that  the  sample  can  be  placed  in  the  instrument,  and  reliably  measured 
without  operator  intervention. 

Several  months  after  receiving  IGAR,  the  LASY-4G  laser  malfunctioned.  As  a 
temporary  measure,  the  laser  was  replaced  with  one  of  Daylight  Solutions’  external 
cavity  quantum  cascade  lasers  (QCL).  This  laser  shifts  the  portion  of  the  spectrum 
which  we  are  using  from  the  LWIR  portion  of  the  spectrum  to  the  MWIR.  The 
laser  currently  being  used  operates  from  5.2  to  5.9pm,  and  can  easily  be  replaced  by 
any  of  the  other  QCLs  in  AFIT’s  possession  [17].  This  allows  the  range  of  useable 
wavelengths  to  be  extended  to  4.84  —  9.7pm.  Furthermore,  this  demonstrates  that 
the  lasers  are  easily  interchangeable.  The  only  attribute  that  limits  this  range  is  the 
responsivity  of  the  pyroelectric  detector.  While  this  specific  pyroelectric  detector  has 
not  been  extensively  tested,  most  detectors  have  a  wide  range  of  useable  wavelengths, 
covering  the  majority  of  the  infrared  portion  of  the  spectrum. 
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Figure  14.  Digram  of  optics  used  in  IGAR. 


Due  to  the  tunable  nature  of  the  instrument,  some  way  to  accurately  and  reliably 
measure  the  wavelength  of  the  laser  was  needed.  For  this,  a  Bristol  721  wavemeter 
[11]  was  used.  The  wavemeter  reports  both  power  and  wavelength  back  to  the  soft¬ 
ware.  Both  of  these  values  are  used  for  the  tuning  of  the  CO2  laser.  The  tuning  of 
the  QCLs  happen  external  to  the  program.  Wavelengths  are  recorded  during  every 
sample  to  ensure  that  the  laser  was  stable  over  the  entire  measurement,  and  the  power 
measurement  is  used  directly  by  the  software  to  calculate  the  DHR  of  the  sample,  as 
detailed  in  Section  6.1. 

Laser  light  is  passed  through  a  series  of  optics  before  illuminating  the  sample, 
as  shown  in  Figure  14.  This  diagram  is  accurate  while  using  the  CO2  laser.  Using 
a  QCL  requires  an  additional  two  mirrors  before  the  beamsplitter  for  alignment.  A 
zinc  selenide  beam  splitter  picks  off  part  of  the  beam  for  spectral  measurements.  A 
50%  duty  cycle  chopper  is  used  in  conjunction  with  a  lock-in  amplifier  to  reduce 
measurement  noise.  A  half-wave  plate  was  included  to  control  the  orientation  of 
the  polarization.  A  replaceable  ND  filter  (0.1,  0.3,  or  1.0)  can  be  placed  to  reduce 
beam  signal  if  necessary.  A  beam  blocker  can  also  be  placed  in  front  of  the  beam 
for  background  measurements  if  necessary.  After  the  filter,  the  light  illuminates  the 
sample.  This  light  is  reflected  and  focused  by  the  mirror. 

Northrop  Grumman  had  Aero  Research  Associates  Inc.  of  Holbrook,  NY,  make 
the  hemi-ellipsoid  mirror  for  IGAR  [7].  This  mirror  has  a  major  axis  of  10.77”  and 
a  minor  axis  of  10.583”.  This  leads  to  a  spacing  of  2”  between  the  upper  and  lower 
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Figure  15.  Diagram  of  the  detector.  The  thick  arrow  is  laser  light’s  vector  reflected 
from  the  sample  off  the  mirror  at  a  high  angle  of  incidence. 

foci,  where  the  detector  and  sample  are  mounted  respectively.  A  slot  0.40”  wide  in 
the  front  of  the  mirror  allows  the  laser  light  into  the  hemisphere  to  illuminate  the 
sample.  Because  the  sample  is  mounted  at  one  focal  point  of  the  mirror,  all  of  the 
light  reflected  by  the  sample  is  reflected  off  the  mirror  to  the  second  focal  point,  where 
the  detector  is  mounted.  The  mirror  has  >  96%  reflectance  throughout  the  1R,  as 
measured  by  NGC  [6],  ensuring  that  nearly  all  of  the  reflected  light  is  collected  and 
delivered  to  the  detector. 

Finally,  the  light  is  collected  by  the  detector.  Northrop  Grumman  designed  a  five¬ 
sided  pyroelectic  detector  for  the  sensor  of  this  system.  The  detector  can  be  seen  in 
Figure  15.  The  power  incident  on  each  face  is  measured  independently  and  reported 
to  the  software.  An  external  controller  is  used  with  the  detector,  and  allows  the 
software  to  specify  a  time  constant  for  the  measurement.  Northrop  Grumman  claims 
that  there  is  some  noise  inherent  to  the  system  due  to  crosstalk  between  the  detector 
channels;  we  at  AFIT  have  not  seen  this  become  an  issue  during  measurements.  This 
effect  is  particularly  apparent  with  specular  samples  when  the  signal  shifts  between 
faces.  Northrop  has  plans  to  enhance  the  detector  by  allowing  each  face  to  be  powered 
and  measured  separately,  eliminating  this  problem. 
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Top  View 


6.3  Alignment 

As  with  any  optical  system,  it  is  crucial  to  have  good  alignment.  IGAR  is  no 
exception.  By  using  the  small  angle  approximation,  any  misalignment  in  the  system 
can  be  described  as  a  combination  of  four  independent  misalignments,  each  with  their 
own  unique  characteristics.  These  four  misalignments  are  x  angular,  x  positional,  y 
angular,  and  y  positional.  A  diagram  of  these  four  cases  can  be  seen  below  in  Figures 
16,  17,  and  18.  Some  of  these  will  adversely  effect  the  performance  of  the  system, 
while  others  can  be  ignored  entirely. 

x  angular  misalignment  is  when  the  beam  is  oriented  at  some  angle  away  from 
the  sample  plane’s  normal,  but  still  parallel  to  the  plane  of  the  optical  bench.  This 
misalignment  will  still  allow  the  equipment  to  capture  uncorrupted  measurements 
from  the  system,  but  there  will  be  an  angular  offset  present  in  the  data.  Due  to 
the  novel  nature  of  this  instrument,  it  is  difficult  to  use  the  data  collected  by  other 
instruments  to  identify  this  misalignment.  In  order  to  identify  this  misalignment,  a 
sample  with  a  distinctive,  angle  dependent  DHR  would  be  needed.  For  instance,  a 
sample  with  a  peak  at  a  known  angle.  Without  this  sort  of  sample,  this  misalignment 
cannot  be  identified.  However,  the  small  angle  approximation  made  earlier  makes 
this  misalignment  negligible. 

The  second  angular  misalignment  is  the  y  angular  misalignment,  as  shown  in 
Figure  17.  This  misalignment  is  negligible  for  small  angles.  For  instance,  consider 
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Top  View 


Top  View 


the  case  where  there  is  a  flat  mirror  in  the  sample  slot.  The  beam  misalignment  will 
simply  move  the  spot  on  the  hemicllipsoidal  mirror  slightly  up  or  down.  Due  to  the 
shape  of  the  mirror,  this  beam  will  still  travel  through  the  second  foci.  For  small 
angles,  this  misalignment  has  no  effect  on  measurements. 

The  other  two  possible  misalignments  are  x  and  y  positional,  as  seen  in  Figures 
16  and  18  respectively.  Both  of  these  are  best  understood  by,  again,  examining 
the  detector  in  image  space,  which  is  the  sample  plane.  The  mirror  does  introduce 
some  aberrations  in  the  system,  which  can  be  ignored  thanks  to  the  small  angle 
approximation.  Positional  y  misalignment  is  the  easiest  to  analyze.  Assuming  again 
that  there  is  only  a  small  misalignment,  the  system  will  be  only  slightly  impacted. 
The  beam  spot  on  the  detector  simply  moves  up  the  face  of  the  detector. 

x  positional  misalignment  is  the  most  difficult  to  analyze  due  to  the  geometry 
of  the  problem.  In  all  of  the  other  misalignments,  the  rotation  of  the  sample  plane 
throughout  the  measurement  can  be  ignored.  Here,  this  effect  cannot  be  ignored; 
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Figure  19.  x  positional  misalignment  at  normal  incidence  and  at  70°. 


rotating  the  sample  plane  means  that  the  x  positional  misalignment  grows.  If  Xq 
is  the  distance  from  the  sample  center  to  the  beam  center  at  normal  incidence,  the 
distance  at  some  angle  6  is  xq/  008(6*).  IGAR  begins  to  make  measurements  at  30°, 
meaning  that  this  length  is  already  15%  larger  when  you  begin  the  measurement.  At 
the  highest  angle,  this  distance  has  grown  to  11  times  the  original  distance.  This  error 
at  both  normal  incidence  and  high  angles  of  incidence  can  be  seen  below  in  Figure 
19. 

If  we  are  analyzing  the  system  with  rays,  this  offset  is  mostly  insignificant.  As 
long  as  the  ray  still  hits  the  detector,  theoretically  we  can  still  measure  the  power. 
With  a  detector  size  of  about  20mm,  we  can  be  off  by  as  much  as  1mm  and  still  make 
measurements.  However,  this  analysis  ignores  the  width  of  the  beam.  If  the  beam  is 
2.5mm  wide,  at  85  degrees,  the  power  of  the  beam  is  distributed  over  28mm.  Even 
with  no  misalignment,  over  20%  of  the  power  is  cast  in  such  a  way  that  it  cannot  be 
seen  by  the  detector.  With  a  misalignment  as  little  as  1mm  at  normal  incidence,  the 
power  lost  increases  to  35%.  This  amount  of  misalignment  can  drastically  impact  the 
accuracy  of  measurements  made  at  high  angles  of  incidence. 
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An  easier  method  for  proper  alignment  might  have  been  including  an  alignment 
mirror.  This  mirror  would  be  oriented  at  45degrees  from  normal  in  the  same  holders 
as  the  calibration  mirrors.  A  retroreflection  would  allow  the  user  to  make  sure  that 
the  laser  is  properly  aligned  to  the  sample  plane,  allowing  much  greater  confidence  in 
the  quality  of  the  measurement. 


6.4  Calibration 


As  shown  above  in  section  6.4,  IGAR  can  be  calibrated  using  the  equation 
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In  this  equation,  pc(0c),  Pwc,j  and  Vdc,j(Ocj )  are  all  known  or  measurable  quantities. 
This  allows  the  equipment  to  easily  calculate  the  value  of  the  calibration  coefficient. 
Actually  making  these  measurements  is  a  bit  more  difficult. 

The  first  term,  pc,j(9c,j )>  is  the  easiest  to  find.  All  that  is  needed  is  the  reflectivity 
of  whatever  sample  we  use  to  calibrate  the  system.  In  practice,  a  gold  coated  mirror 
is  used  for  calibration.  This  mirror  has  high  reflectivity  over  the  entire  IR  portion 
of  the  spectrum,  about  97%.  It  is  important  to  note  that  the  angle  of  incidence 
is  included  in  this  term.  This  has  less  to  do  with  the  reflectivity  changing  as  a 
function  of  angle  as  it  does  with  the  orientation  of  the  mirror  during  calibration.  The 
angular  dependence  of  this  sample  can  be  calculated  using  Frensel’s  equations.  For 
other  samples,  a  measurement  of  the  angular  characteristics  of  the  calibration  sample 
would  need  to  be  known. 

The  next  term  is  Pwc,j  ■  This  is  simply  a  measurement  of  the  incident  power  during 
the  calibration.  The  wavemeter  can  easily  make  this  measurement  and  relate  it  to 
the  system  for  calculations.  The  incident  power  is  measured  as  a  function  of  face  to 


50 


allow  for  the  power  of  the  laser  to  drift  even  during  calibration.  This  allows  an  even 
more  robust  system  with  no  cost  to  either  complexity  or  measurement  time. 

Finally,  there  is  the  term  Vdc,j{Oc,j)-  There  are  three  important  aspects  of  this 
measurement:  the  voltage,  the  angle  of  incidence  and  the  face  being  measured.  This 
calibration  method  allows  each  face  to  be  calibrated  individually,  but  comes  with  the 
cost  that  all  of  the  reflected  power  is  incident  onto  a  single  face.  Without  this  being 
true,  the  reflectivity  that  the  detector  sees  would  be  distributed  between  multiple 
faces,  meaning  that  the  pc  term  is  not  actually  known.  By  controlling  the  orientation 
of  the  sample,  the  system  can  be  adjusted  such  that  all  of  the  light  incident  on  the 
calibration  sample  will  be  reflected  to  the  face  of  the  detector  being  calibrated. 

When  this  assumption  holds  true,  the  power  incident  on  the  calibration  mirror, 
the  reflectivity  of  the  mirror,  and  the  voltage  of  the  detector  for  the  incident  power 
are  all  known  quantities.  With  these  measurements,  the  system  can  be  calibrated, 
and  used  for  making  DHR  measurements. 

Originally,  Northrop  Grumman  designed  IGAR  to  be  calibrated  before  and  after 
making  a  set  of  measurements.  This  was  clue  to  the  fact  that  initially,  IGAR  had  no 
way  to  measure  the  power  of  the  laser  during  the  measurement.  Originally,  IGAR 
was  designed  to  omit  the  Pw  and  Pwc  terms  from  the  DHR  equation  and  calibration 
equation.  Instead,  the  power  was  measured  before  the  initial  calibration  and  after 
the  final  calibration.  If  these  two  numbers  were  the  same,  the  power  could  simply  be 
ignored.  Pw  =  Pwc,  and  could  be  eliminated  from  the  equation.  If  these  two  numbers 
were  close,  a  linear  gradient  was  applied  to  the  data.  This  effectively  functioned  as 
a  Pw  that  increased  from  the  initial  value  to  the  final  value  over  the  course  of  the 
measurement.  By  adding  the  wavemeter  to  the  system,  the  incident  power  could  be 
measured  in  real  time,  which  eliminated  the  need  for  this  scaling  effect,  and  allowed 
the  final  calibration  to  be  ignored. 
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Eventually,  further  testing  revelaved  the  initial  calibration  could  be  forgone  as 
well.  Calculating  the  calibration  coefficient  as  done  in  Equation  (86)  gave  the  system 
the  flexibility  needed  to  do  this.  By  including  the  power  term,  any  short  or  long  term 
variations  in  the  laser  power  could  be  accounted  for.  The  responsivity  of  the  detector 
does  not  change  (except  as  a  function  of  wavelength),  which  makes  this  term  more  or 
less  a  constant.  Now,  calibration  is  only  performed  as  a  verification  method,  or  when 
the  wavelength  of  the  laser  is  changed. 

Note  that  there  is  some  small  wavelength  dependence  in  the  responsivity  of  the  de¬ 
tector.  A  pyroelectric  detector  functions  on  the  change  in  temperature  of  the  material 
covering  the  detector,  making  the  detector  responsivity  a  function  of  the  absorptance 
of  the  material  used  to  cover  the  detector.  In  practice,  the  material  covering  the 
detector  is  a  dark,  very  absorptive  material,  making  this  dependence  on  wavelength 
negligible.  Even  without  this  assumption,  this  responsivity  term  is  absorbed  into  the 
calibration  coefficient,  and  can  be  calibrated  out  of  the  measurement. 

At  one  point,  Northrop  Grumman  considered  using  a  diffuse  sample  to  calibrate 
all  five  detector  faces  simultaneously.  The  idea  was  that  a  diffuse  sample  would 
illuminate  each  face  identically.  Knowing  the  incident  power,  and  that  this  incident 
power  was  equally  distributed  between  all  five  faces,  the  calibration  constants  could 
be  calculated  simultaneously.  In  theory,  this  could  work:  knowing  the  reflectivity 
of  a  diffuse  sample  and  the  size  of  the  faces,  the  incident  power  could  easily  be 
calculated.  However,  there  are  no  well  defined  diffuse  DHRs  in  the  LWIR  portion  of 
the  spectrum.  At  shorter  wavelengths,  a  sample  like  Spectralon  or  Infragold  would 
function  well.  However,  there  are  no  diffuse  NIST  traceable  standards  around  the 
10pm  wavelength  that  could  be  used  to  calibrate  the  machine.  Northrop  instead 
opted  for  the  multiple  mirror  calibration  methodology  previously  described. 
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6.5  Software 


All  of  the  software  for  IGAR  was  written  in  Labview.  The  majority  of  the  code 
is  used  to  facilitate  the  use  of  the  device,  controlling  parameters  such  as  angle  of 
incidence  and  operating  the  pyroelectric  detector.  However,  there  are  a  few  sections 
of  the  code  which  are  worth  discussing  here.  One  of  these  is  the  laser  autotuning 
subroutine.  This  subroutine  is  used  for  automatically  tuning  the  CO2  laser  to  the 
desired  wavelength  and  optimizing  the  power.  Creating  an  algorithm  to  do  this 
turned  out  to  be  nontrivial  due  to  nonlinearities  in  the  laser.  For  instance,  the 
optimal  position  of  the  micrometer  used  for  tuning  the  laser  would  shift  depending 
on  the  temperature  of  the  laser  housing.  Another  significant  portion  of  the  code 
is  used  to  calculate  the  power  incident  on  the  pyroelectric  detector.  This  algorithm 
analyzes  the  data  on  all  five  faces  to  determine  which  ones  will  yield  the  most  accurate 
measurement. 
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VII.  Algorithms 


There  are  two  notable  algorithms  implemented  in  this  software  that  are  used  by 
IGAR  to  enable  and  aid  measurements.  The  first  of  these  is  the  Laser  Autotuning 
algorithm.  This  algorithm  is  an  automated  sequence  which  is  used  by  the  software 
to  control  the  CO2  laser  previously  mentioned.  This  method  allows  the  equipment 
to  tune  the  laser  to  the  desired  wavelength  without  operator  intervention,  allowing 
the  system  to  make  DffR  measurements  of  a  single  sample  at  multiple  wavelengths 
efficiently.  The  second  is  the  algorithm  mentioned  in  Section  5.7.  This  algorithm 
intelligently  selects  which  faces  are  used  to  calculate  the  DHR  in  order  to  reduce  the 
uncertainty  of  the  measurement.  Both  of  these  algorithms  were  written  in  Labview 
and  have  been  extensively  tested  in  laboratory  work. 

7.1  Laser  Autotuning  Algorithm 

The  laser  auto  tuning  subroutine  was  designed  to  automatically  tune  the  CO2 
laser  to  the  desired  wavelength  without  operator  intervention.  The  CO2  laser  is 
a  grating-tuned  laser  with  a  computer  controlled  micrometer  moving  the  grating. 
Tuning  the  laser  to  the  desired  wavelength  proved  to  be  nontrivial.  The  peak  power 
and  exact  position  of  the  micrometer  for  a  single  wavelength  tended  to  shift  depending 
on  external  conditions.  The  exact  position  varied,  but  by  a  small  enough  amount  that 
the  approximate  position  of  the  wavelength  is  known. 

Knowing  this  approximate  position,  1  designed  an  algorithm  to  tune  the  laser  to 
the  peak  power  at  the  desired  wavelength.  The  diagram  of  this  process  can  be  seen 
in  20.  The  subroutine  functions  as  follows.  First,  the  micrometer  is  moved  just  short 
of  the  approximate  position  of  the  peak  power.  The  micrometer  is  then  incremented 
forward  until  lasing  is  observed.  At  this  point,  the  wavelength  is  checked  by  the 
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wavemeter.  Assuming  the  wavelength  is  the  one  we  are  looking  for,  the  subroutine 
continues  moving  the  micrometer  forward  until  the  power  of  the  laser  begins  to  drop. 
At  this  point,  the  micrometer  backs  up,  and  begins  moving  forward  again  with  smaller 
steps.  This  process  is  repeated  three  times.  Once  the  peak  is  identified,  the  microm¬ 
eter  resets,  and  moves  to  the  recently  measured  peak  position.  The  wavelength  is 
checked  once  again  and  recorded.  The  entire  process  takes  about  five  minutes,  and 
has  been  extensively  tested.  It  will  reliably  find  the  peak  power  with  no  operator 
intervention,  allowing  us  to  run  a  single  sample  at  multiple  wavelengths. 


|  SI  |  Set  initial  microme¬ 
ter  position 

Initialize  micrometer,  move  to  predefined  initial  position. 

|  S2  |  Search  for  lasing 

Increment  micrometer  until  lasing  is  measured.  If  the  laser  is  not 
lasing,  end  subroutine  and  notify  user. 

[~S3~|  Check  A 

Measure  wavelength  via  Bristol  721  Wavemeter,  and  compare  with 
predefined  values.  If  the  wavelength  is  not  correct,  end  subroutine 
and  notify  user. 

S4  Pause 

Wait  for  laser  to  stabilize. 

S5  Power  loop 

Increment  micrometer  position  until  peak  power  is  observed.  After 
peak,  decrement  the  micrometer  position,  and  re-find  peak  power 
using  a  smaller  step  size.  This  is  repeated  with  three  step  sizes.  If 
the  peak  cannot  be  found,  notify  the  user  and  end  the  subroutine. 

S6  Pause 

Wait  for  laser  to  stabilize. 

|  ST  |  Recheck  power 

Reset  the  micrometer,  move  to  peak  power  location  and  remeasure 
power.  If  the  power  is  not  the  same  as  the  power  measured  in  S5, 
return  to  SI. 

S8  Complete 

Report  power  and  wavelength  to  software. 

Figure  20.  Algorithm  for  automatically  tuning  the  CO2  laser  to  the  desired  wavelength 
without  operator  intervention. 
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7.2  IGAR  Algorithm 


The  other  notable  piece  of  code  is  the  algorithm  for  calculating  the  power  reflected 
off  the  sample  using  the  measurements  from  the  pyroelectric  detector.  To  capture  the 
data,  the  software  takes  measurements  at  a  variety  of  time  constants  with  a  variable 
number  of  measurements.  The  mean  voltage  and  standard  deviation  are  calculated, 
which  are  converted  into  power  measurements  using  the  calibration  data,  as  detailed 
in  Chapter  V.  The  algorithm  is  then  called  to  calculate  the  most  accurate  power 
measurement  using  the  data  from  all  five  faces. 

The  algorithm,  designed  by  Northrop  Grumman,  works  based  on  the  assumption 
that  not  all  faces  are  contributing  data  to  the  measurement.  For  instance,  if  all  of 
the  power  is  incident  only  on  one  face,  four  of  the  faces  are  contributing  nothing 
but  noise  to  the  measurement.  By  carefully  ignoring  these  noncontributing  faces,  the 
measurements  can  be  made  more  accurate.  The  algorithm  operates  in  this  manner, 
and  can  be  seen  in  Figure  21. 

To  begin,  the  algorithm  calculates  the  power  and  uncertainty  using  all  five  faces. 
Both  of  these  are  simply  the  sums  of  their  respective  parts.  Next,  the  algorithm 
drops  the  face  with  the  lowest  power  from  set.  Power  and  uncertainty  are  recalculated 
using  this  reduced  set.  These  values  are  then  compared  to  the  previous  power  and 
uncertainty  measurements.  If  the  new  power  is  within  one  uncertainty  of  the  previous 
power,  it  is  assumed  that  the  dropped  face  contributed  nothing  to  the  measurement. 
The  process  is  then  repeated  by  dropping  another  face,  until  the  power  is  outside  the 
uncertainty.  The  last  power  within  the  uncertainty  is  reported  back  to  the  program 
and  used  to  calculate  the  DHR. 

Both  of  these  algorithms  have  been  successfully  implemented  in  IGAR  and  used 
to  make  measurements.  The  autotuning  algorithm  has  been  implemented  and  used 
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Figure  21.  Algorithm  for  determining  most  accurate  DHR  from  five  sided  detector. 

successfully  to  make  DHR  measurements  on  a  single  sample  at  multiple  wavelengths. 
For  the  samples  examined  in  this  paper,  only  a  single  wavelength  was  used,  but  IGAR 
has  the  capability  of  making  multiple  measurements  automatically.  The  IGAR  algo¬ 
rithm  has  been  successfully  used  on  all  of  the  measurements  shown  in  the  next  section. 
The  outputs  of  this  algorithm  have  been  compared  to  the  results  without  using  the 
algorithm,  and  been  shown  that  the  algorithm  decreases  measurement  uncertainty 
while  retaining  the  accuracy  of  the  DHR  given  the  measured  data. 


57 


VIII.  Measurements 


IGAR  has  taken  a  number  of  measurements  of  different  samples.  Two  of  important 
note  are  an  aluminum  block  with  significant  surface  roughness  and  a  measurement  of 
the  gold  calibration  mirror.  The  aluminum  block  is  a  sample  with  high  reflectivity 
which  is  very  diffuse.  This  means  that  the  light  reflected  by  the  sample  will  be 
scattered  over  all  five  faces  of  the  detector.  The  gold  mirror  is  a  much  more  specular 
sample,  which  will  show  the  reflection  walking  over  the  detector. 

For  the  following  measurements,  one  of  the  QCLs  was  used  to  illuminate  the 
sample.  The  laser  was  tuned  to  5.504pm.  This  particular  laser  emits  vertically 
polarized  light  (parallel  to  the  y  axis),  and  no  optics  which  affect  polarization  were 
used.  Both  of  these  measurements  were  taken  using  the  same  calibration  data,  which 
was  taken  slightly  before  the  measurement. 

8.1  Aluminum  Block 

The  aluminum  block  is  best  compared  to  the  Black  Oxidized  Steel  scatterer  shown 
in  Section  4.5.  While  the  samples  have  unequal  DHR  measurements,  they  both  share 
a  number  of  characteristics.  Before  continuing,  we  must  address  the  peak  at  the  end 
of  the  DHR  measurement,  seen  in  Figure  22.  This  peak  occurs  due  to  a  misalignment 
of  the  incident  laser  beam.  The  fourth  face  of  the  detector  could  not  be  properly 
calibrated,  which  makes  the  contribution  of  this  face  to  the  DHR  greater  than  it 
should  be.  Ignoring  this  error,  the  voltage  seen  by  each  face  of  the  detector  is  very 
similar  to  what  would  be  expected  given  a  diffuse  scatterer. 

This  is  shown  in  Figure  23.  While  the  Black  Oxidized  Steel  scatterer  is  not 
completely  accurate  to  the  aluminum  block,  the  signals  across  each  face  are  consistent 
with  each  other.  Face  2  seems  to  have  the  greatest  discrepancy,  but  this  is  due  to  the 
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Aluminum  Block  2 


Figure  22.  The  DHR  and  face  voltages  for  the  rough  surfaced  aluminum  block. 


normalization.  Examining  Figure  22,  this  face  has  the  lowest  voltage  of  the  five.  The 
discrepancy  in  Figure  23  is  due  to  poor  normalization.  The  other  four  faces  follow 
similar  patterns  to  the  Black  Oxidized  Steel  scatterer.  Faces  1  and  3  are  the  Left  and 
Right  faces,  respectively,  so  it  makes  sense  that  the  signal  gradually  drops  off  these 
faces  as  the  sample  approaches  grazing  angles.  When  this  happens,  the  signal  shift 
to  the  Front  face  (Face  4)  of  the  detector,  which  happens  in  both  measurements  and 
the  simulation. 

With  some  further  simulations,  the  aluminum  scatterer  could  be  more  accurately 
modeled,  and  simulation  could  be  made  more  accurate  to  what  was  measured  in  lab. 
Nonetheless,  the  correlation  between  the  simulation  of  a  diffuse  scatterer  and  the 
measured  data  from  IGAR  shows  that  both  systems  work  well  in  measuring  DHR. 

8.2  Gold  Mirror 

At  first  glance,  the  face  measurements  in  Figure  25  do  not  appear  consistent  with 
one  another.  This  is  again,  due  to  normalization.  Faces  1,  2  and  3  have  very  little 
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Figure  23.  Simulated  measurements  using  the  Black  Oxidized  Steel  scatterer.  Dashed 
red  line  is  the  voltage  measurements  of  the  rough  aluminum  block. 

signal  compared  to  Faces  4  and  5.  Because  both  the  gold  mirror  and  Gold  Metallic 
scatterer  are  such  specular  samples,  nearly  all  of  the  energy  is  concentrated  on  the 
front  and  top  faces  of  the  detector.  The  actual  voltages  can  be  seen  in  Figure  24. 
The  numbers  from  the  simulation  are  similarly  low. 

Faces  4  and  5  hold  the  most  information  in  this  sample.  What  is  seen  in  Figure  25 
is  that  the  power  reflected  by  the  sample  transition  from  the  top  face  to  the  front  face 
of  the  detector.  This  is  seen  happening  gradually  in  the  case  of  the  Gold  Metallic 
scatterer,  and  very  quickly  in  the  case  of  the  gold  mirror.  This  is  mainly  due  to 
the  construction  of  the  simulation.  In  the  simulation,  the  sample  was  square,  and 
equally  illuminated  over  its  entire  surface.  For  the  measurement,  the  gold  mirror  was 
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Mirror  4 


Figure  24.  The  DHR  and  face  voltages  for  the  gold  calibration  mirror. 


illuminated  by  a  Gaussian  laser  beam.  This  led  to  more  power  being  concentrated  at 
the  middle  of  the  beam,  which  explains  the  sudden  dropoff  between  faces. 

8.3  Experimental  Results 

These  results  show  that  the  underlying  theory  of  this  experiment  is  accurate.  The 
theoretical  five-sided  detector  gives  results  that  are  similar  to  what  can  be  measured 
in  lab.  The  diffuse  black  oxidized  steel  scatterer  gave  results  that  very  accurately 
resembled  the  aluminum  block  measured  by  IGAR.  With  a  more  accurate  model,  the 
simulated  output  of  the  system  would  be  even  closer  to  real  results.  Likewise,  the 
gold  scatterer  followed  the  same  patterns  found  in  the  gold  mirror.  Because  both 
scatterers  were  highly  specular,  their  results  followed  more  closely,  regardless  of  the 
inaccuracies  of  the  underlying  model.  With  some  further  time  spent  building  models 
from  measured  BRDF,  I  am  confident  that  the  simulation  of  this  system  could  even 
more  closely  resemble  the  experimental  data. 
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Figure  25.  Simulated  measurements  using  the  Gold  Metallic  scatterer.  Dashed  red  line 
is  the  voltage  measurements  of  the  gold  calibration  mirror. 
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IX.  Conclusions  and  Future  Work 


9.1  Conclusions 

As  our  results  have  shown,  IGAR  can  make  measurements  of  the  DHR  of  a  sample. 
This  ability  stems  from  the  design  of  the  five  sided  pyroelectric  detector.  While  the 
detector  cannot  capture  all  of  the  light  reflected  from  the  sample,  the  error  of  the 
measurement  falls  within  the  uncertainty  of  the  sample,  making  this  source  of  error 
negligible. 

A  robust  calibration  methodology  ensures  that  the  system  makes  accurate  mea¬ 
surements.  By  carefully  designing  the  system,  the  system  can  easily  and  accurately 
calibrated.  This  calibration  methodology  has  the  added  benefit  of  being  able  to  be 
calibrated  only  knowing  the  exact  characteristics  of  the  calibration  mirrors.  All  other 
elements  of  the  system,  such  as  the  ratio  of  the  beamsplitter  or  the  voltage  to  power 
ratio  of  the  detector  faces  can  be  calibrated  out.  This  allows  for  high  accuracy  and 
ease  of  use. 

Many  of  the  errors  found  in  the  system  can  be  traced  back  to  the  misalignment 
of  the  laser  beam.  With  a  better  alignment  procedure,  many  of  these  errors  can  be 
eliminated.  Even  so,  the  measurements  taken  by  the  system  are  promising.  Many 
of  the  attributes  found  in  the  simulation  of  different  scatterers  can  be  backed  up 
by  experimental  data.  With  more  accurate  scattering  models,  the  experimentally 
collected  data  can  be  further  verified. 

IGAR  shows  much  promise  as  a  system  capable  of  making  DHR  measurements  at 
high  angles  of  incidence.  With  it’s  novel  detector,  tunable  laser  capabilities,  and  ro¬ 
bust  calibration  procedures,  it  can  make  accurate  measurements  of  DHR  and  increase 
the  already  impressive  measurement  capabilities  of  the  United  States  Air  Force. 
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Appendix  A.  Measurement  Uncertainty  Expanded 


In  this  appendix,  we  go  from  RD  eqn  whatever  to  RD  6.21.  Eqn  6.21.  Let’s  begin 
with 

(87) 

3= 1 

Assuming  the  system  is  perfectly  calibrated  (i.e.,  there  is  no  uncertainty  in  the  calibra¬ 
tion  coefficient),  there  are  three  terms  with  uncertainties:  uncertainty  in  the  voltage 
measured  at  the  detector  (A Vdj),  power  measured  by  the  wavemeter  (A Pw),  and  the 
uncertainty  of  the  measurement  of  the  DHR  of  the  sample  (A p(0)).  Accounting  for 
these  uncertainties,  the  equation  above  becomes 


Pifi)  +  Ap(0) 


r  VdM  +  &Vdj(0) 

4?  c’j  PW  +  APW 

3= 1 


(88) 


Expanding  the  voltage  term, 


5 

m + a  p(e)  =  y  ccj 

l=i 


P  +  A  P 

1  W  I  ‘ — 1 1  w 


c’jPw  +  APw 

3  =  1 


(89) 


Expanding  the  power  term, 


m + Ap(#>  =  E  g, 

l=i 


VaM 


•3  p 

1  II 


+  Cc’- 


1=1 


Y*M 

pi 


A  P, 


W 

5  a  t  r  /  r\\  5 


1=i 


AVdJ(6) 


c>3  p 

'll. 


y  ^c.i 

1=1 


AhdiJ(0) 


P2 

W 


A/> 


(90) 


The  uncertainty  cross  terms  can  be  ignored.  The  term  without  any  uncertainty  is 
what  was  previously  defined  as  the  DHR  of  the  measurement,  and  can  be  used  to 
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cancel  out  the  p{9)  term.  This  leaves 


5 


Z—/  C,J  P2 
3= 1  W 


c’3  p 
3= 1  U 


(91) 


Using  the  definition  of  DHR,  the  equation  can  be  further  simplified  down  to 


A  P  _  , 

Ap(0)  =  M#)^  +  Ec<=, 


p„ 


3= 1 


AM) 

p 

7/? 


(92) 


At  this  point,  we  see  that  even  allowing  uncertainty  in  calibration  will  not  effect  the 
uncertainty  of  the  entire  measurement.  This  uncertainty  term  is  multiplied  by  the 
uncertainty  of  the  voltage  measurement,  and  can  therefore  be  ignored. 

Note  that  the  first  term  is  independent  of  face,  while  the  second  term  is  not. 
Because  of  this,  ignoring  the  contributions  of  some  of  the  faces  can  decrease  the 
uncertainty  of  the  measurement.  This  is  the  motivation  for  section  whatever  where 
we  talk  about  the  1GAR  algorithm. 
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Appendix  B.  MATLAB  Code 


Below  is  the  MATLAB  code  used  for  the  simulations  found  above  in  Chapter 
IV.  <Lact  is  found  by  running  code  “cookTorrance”.  <f>det  is  found  running  the  code 
“detectorOptimal3” .  Plots  covering  a  range  of  angles  were  calculated  by  running 
this  function  over  this  angular  range  and  plotting  the  results.  “cookTorrance VL” 
functions  the  same  as  “cookTorrance” ,  but  has  input  vector,  as  opposed  to  the  angle 
of  the  incident  light,  and  is  called  by  “detectorOptimal3” . 

2.1  cookTorrance. m 

function  [R,DHR]  =  cookTorrance (P,T) 

if  nargin  <  2 
P  =  45; 

T  =  45; 

end 


7.  [sample  name] 

> 

FO 

m 

7.  acrylic  blue 

> 

0.035 

0.037 

7.  black  oxidized  steel 

> 

0.035 

0.19 

7.  black  plastic  soft 

> 

0.136 

0.325 

7.  blue  rubber 

> 

0.0366 

0.276 

7.  ball  gold  metllic  2 

> 

0.629 

0.144 

7.  ipswitch  pine 

> 

0.0192 

0.12 

7.  metallic  blue 

> 

0.638 

0.231 

7.  metallic  silver 

> 

0.661 

0.195 

7.  nickel 

> 

0.733 

0.0424 

7.  provincial 

> 

0.114 

0.0864 

7.  red  oak 

> 

0.0248 

0.116 

FO  =  0.629;  m  =  0.144 

n  =  (1+sqrt (FO) ) /(1-sqrt (FO) ) ; 


p  =  512; 

[phi  theta]  =  meshgrid(0:pi/2/p:pi/2,  0 : 2*pi/p : 2*pi) ; 
7.7.  geometry  calculations 
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[x  y  z]  =  sph2cart (theta, phi , ones (size (theta) )) ; 

V(l,  : ,  =  x; 

V (2, :  ,  =  y; 

V (3 , :  ,  =  z; 

L  =  [cosd(P) *cosd(T)  sind(P) *cosd(T)  sind(T)] ; 

Lp ( 1 ,  :  ,  : )  =  L(l)*ones (size(V(l , :,:))); 

Lp (2 , : , : )  =  L (2) *ones (size (V(l ,:,:))) ; 

Lp (3 , : , : )  =  L (3) *ones (size (V(l ,:,:))) ; 

temp  =  V  +  Lp; 

HI  =  sqrt (temp(l ~2  +  temp (2, : , : ) . ~2  +  temp (3, : , : ) . ~2) ; 

H(1 ,  :  ,  : )  =  temp(l , : , ./HI; 

H(2, : , : )  =  temp(2, : , : ) . /HI ; 

H (3 , : , : )  =  temp (3 , : , : ) . /HI ; 

cosTheta  =  V(1 *H(1 )  +  V(2 *H(2 )  +  V(3, : , : ) . *H(3, : , : ) 

•/.•/.  cook-torrance  constants 

cosDelta  =  H(3, : ,  : ) ; 
sinDelta  =  sqrt(l-H(3, : , : ) . ~2) ; 
tanDelta  =  sinDelta. /cosDelta; 

D  =  l/m~2/cosDelta. ~4.*exp(-(tanDelta/m) . ~2) ; 

Gp ( 1 , : , : )  =  ones (size (V(l ,:,:))) ; 

Gp(2,:,:)  =  2*H(3 *V(3, /cosTheta; 

Gp(3,:,:)  =  2*H(3 ,:,:). *Lp(3 ,:,:). /cosTheta; 

G  =  min(Gp,  []  ,  1)  ; 

•/.•/.  fresnel  calculations 

n  =  10.39993  +  i*40; 
g  =  sqrt(n~2  +  cosTheta. ~2  -1); 

F  =  1/2* (g-cosTheta) . ~2 . / (g+cosTheta) . ~2 

.*(1  +  (cosTheta. * (g+cosTheta)  -  1).~2 
./ (cosTheta. * (g-cosTheta)  +  1).~2); 

•/.'/.  cook-torrance  calculations 

K  =  l/pi*D . *G . /Lp(3 /V (3 *F; 

if  1 
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°/„[x  y  z]  =  sph2cart  (theta, phi ,  ones  (size  (theta) ))  ; 
surf (x,y,z, squeeze (abs(K)) , ’ EdgeAlpha’ ,0) 

%surf (x , y , z , z , ’ EdgeAlpha ’ , 0) 
axis  equal 

axis ( [-1  1-11-0  1] *1.2) 

end 

nanFinder  =  isnan(K) ; 

S  =  size (K) ; 
for  j  =  1 : S(l) 

for  k  =  1 : S (2) 

if  nanFinder(j ,k)  ==  1 
K(j,k)  =  0; 

end 

end 

end 

R  =  squeeze (K) ; 

DHR  =  sum(sum(R*pi/p*pi/p*4 . *sin(pi/2-phi) . *cos (pi/2-phi) ) ) 

if  nargin  <  2 

R  =  0;  DHR  =  0; 

end 

end 

2.2  cookTorranceVL.m 

function  [R]  =  cookTorranceVL(V,Lp) 


y„  [sample  name] 

°/„  black  plastic  soft 
%  blue  rubber 
%  ball  gold  metllic  2 
%  ipswitch  pine 
%  metallic  blue 
°/„  metallic  silver 
°/„  nickel 
%  provincial 
°/„  red  oak 


> 

F0 

m 

> 

0.136 

0.325 

> 

0.0366 

0.276 

> 

0.629 

0.144 

> 

0.0192 

0.12 

> 

0.638 

0.231 

> 

0.661 

0.195 

> 

0.733 

0.0424 

> 

0.114 

0.0864 

> 

0.0248 

0.116 

F0  =  0.661;  m  =  0.195; 


n  =  (1+sqrt (FO) ) /(1-sqrt (FO) ) ; 


%%  geometry  calculations 

VI  =  sqrt (V(l ~2  +  V(2,:,:).~2  +  V(3, : , : ) . ~2) ; 

V(l,  :  ,  :)  =  V(l, : , :) ./VI; 

V (2 ,  :  ,  :)  =  V (2 , : , :) ./VI; 

V (3 , : , :)  =  V (3 , : , :) ./VI; 

LI  =  sqrt (Lp ( 1 ~2  +  Lp(2,:,:).~2  +  Lp(3 “2  ); 

Lp ( 1 , : , :)  =  Lp(l, : , :) ./LI; 

Lp (2 , : , :)  =  Lp(2, : , :) ./LI; 

Lp (3 , : , :)  =  Lp(3, : , :) ./LI; 

temp  =  V  +  Lp; 

HI  =  sqrt (temp(l ~2  +  temp (2, : , : ) . ~2  +  temp (3, : , : ) . ~2) ; 

H(l,  :  ,  :)  =  tempCl, : , :) ./HI; 

H(2, : ,  :)  =  temp (2, : , :) ./HI; 

H (3 , : , : )  =  temp (3 , : , : ) . /HI ; 

cosTheta  =  V(1 *H(1 )  +  V(2, : , : ) . *H(2, : , : )  +  V(3, : , : ) . *H(3, : , : ) ; 
•/.'/.  cook-torrance  constants 
cosDelta  =  H (3 

sinDelta  =  sqrt(l-H(3, : , : ) . ~2) ; 
tanDelta  =  sinDelta. /cosDelta; 

D  =  l/m~2/cosDelta. ~4.*exp(-(tanDelta/m) . ~2) ; 

Gp ( 1 , : , : )  =  ones (size ( V ( 1 ,:,:))) ; 

Gp(2,:,:)  =  2*H(3 *V(3, /cosTheta; 

Gp(3,:,:)  =  2*H(3 *Lp(3 /cosTheta; 

G  =  min(Gp,  []  ,  1) ; 

°/o7o  fresnel  calculations 

g  =  sqrt(n~2  +  cosTheta. "2  -1); 

F  =  1/2* (g-cosTheta) . ~2 . / (g+cosTheta) . "2 

.*(1  +  (cosTheta. * (g+cosTheta)  -  1).~2 
./ (cosTheta. * (g-cosTheta)  +  1).~2); 

U  cook-torrance  calculations 
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K  =  l/pi*D . *G . /Lp(3 ,:,:). /V (3 ,:,:). *F; 


if  0 

[x  y  z]  =  sph2cart (theta, phi , ones (size (theta) )) ; 
surf (x,y,z, squeeze (abs(D)) , ’ EdgeAlpha’ ,0) 
axis  vis3d 

axis ([-2  2-2202]) 

end 

nanFinder  =  isnan(K) ; 

S  =  size (K) ; 
for  j  =  1 : S(l) 

for  k  =  1 : S (2) 

if  nanFinder(j ,k)  ==  1 
K(j,k)  =  0; 

end 

end 

end 

R  =  K;  /(squeeze (F)  ; 

°/„DHR  =  sum(sum(R*pi/p*pi/p*4 .  *sin(pi/2-phi)  .  *cos  (pi/2-phi) ) )  ; 
end 

2.3  detectorOptimal3.m 

function  [S,  PhiVec,  PhiStruct]  =  detector0ptimal3 (theta,  phi) 

if  nargin  <  1 
theta  =  35; 

end 

if  nargin  <  2 
phi  =  0; 

end 


%%  definitions 
nl  =  1;  n2  =  4; 


xEMax  =  1/2;  %  1/2  width 
dxE  =  0.1; 

yEMax  =  1/2;  %  1/2  width 
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dyE  =  0.1; 


xDMax  =  10;  %  1/2  width 
dxD  =  0.1/1; 

yDMax  =  20;  %  1/2  width  <-  longer  dimension 
dyD  =  0.1/1; 

zDMax  =  10;  %  full  height 
dzD  =  0.1/1; 

[xe  ye]  =  meshgrid (-xEMax:dxE:xEMax,  -yEMax:dyE:yEMax) ; 

Lo  =  [cosd(theta) *cosd(phi)  cosd(theta) *sind(phi)  sind(theta)] 
N  =  [0  0  1]; 

%%  top  face 

%  z  held  constant,  integrate  over  x  and  y 
t  =  []  ;  f  =  []  ;  1  =  []  ; 

[t.xd  t.yd]  =  meshgrid (-xDMax :dxD: xDMax,  -yDMax : dyD: yDMax) ; 
[f.yd  f.zd]  =  meshgrid (-yDMax: dyD: yDMax,  0+0 . 0001 :dzD : zDMax) ; 
[l.xd  l.zd]  =  meshgrid (-xDMax :dxD: xDMax,  0+0 . 0001 :dzD : zDMax) ; 

S  =  size (xe) ; 

PhiTop  =  zeros (size (t .xd) ) ; 

PhiFront  =  zeros (size (f .yd)) ; 

PhiLeft  =  zeros(size(l.xd)) ; 

PhiBack  =  zeros (size (f .yd)) ; 

PhiRight  =  zeros (size (1 .xd) ) ; 

if  1 

for  K  =  1 : S(l) 

for  L  =  1 : S (2) 

if  (abs(xe(K,L))  <  xDMax)  &&  (abs(ye(K,L))  <  yDMax) 

Ve  =  [xe(K,L)  ye(K,L)  0] ; 

%°/„  top  face 
clear  V  Vd  Lp 
Vd(l , : , : )  =  t.xd; 

Vd(2, : , : )  =  t.yd; 

Vd(3,:,:)  =  zDMax*ones (size (t . xd) ) ; 

V(l,  :  ,  :)  =  Vd(l , : , :)  -  Ve(l) ; 

V (2 , : , :)  =  Vd(2, : , :)  -  Ve(2) ; 
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V (3 , =  Vd(3, : , :)  -  Ve(3); 


Lp ( 1 ,  :  ,  :  )  =  Lo (1) *ones (size(V(l , :,:))); 

Lp(2,:,:)  =  Lo (2) *ones (size ( V ( 1 ,:,:))) ; 

Lp(3,:,:)  =  Lo (3) *ones (size ( V ( 1 ,:,:))) ; 

R  =  sqrt (V(l ~2  +  V(2,:,:).~2  +  V(3,:,:).~2) 
cosThetaE  =  V (3 /R ; 
cosThetaD  =  V(3,:,:)./R; 

A  =  ref lection(nl ,n2 , cosThetaD) ; 

I  =  cookTorranceVL(V,Lp) ; 

Lr  =  A. *1 . *cosThetaE; 

PhiTop  =  PhiTop  +  squeeze (Lr . *A . *cosThetaD 
. *1 . /R. ~2*dxE*dyE*dxD*dyD) ; 


U  front  face 
clear  V  Vd  Lp 

Vd(l,:,:)  =  xDMax*ones(size(f .yd)) ; 
Vd(2, : , : )  =  f .yd; 

Vd(3, : , : )  =  f .zd; 

V(l, : , :)  =  Vd(l, : , :)  -  Ve(l); 

V (2 , :,:)  =  Vd(2, : ,:)  -  Ve( 2); 

V (3 , :,:)  =  Vd(3 , : , :)  -  Ve(3); 

VFront  =  V; 


Lp ( 1 , : , : )  =  Lo (1) *ones (size(V(l , :,:))); 

Lp(2,:,:)  =  Lo (2) *ones (size ( V ( 1 ,:,:))) ; 

Lp(3,:,:)  =  Lo (3) *ones (size ( V ( 1 ,:,:))) ; 

R  =  sqrt (V(l ,:,:). ~2  +  V(2,:,:).~2  +  V(3,:,:).“2) 
cosThetaE  =  V(3,:,:)./R; 
cosThetaD  =  V(l,:,:)./R; 

A  =  ref lection(nl ,n2 , cosThetaD) ; 

I  =  cookTorranceVL(V,Lp) ; 

Lr  =  A. *1 . ^cosThetaE; 

PhiFront  =  PhiFront  +  squeeze (Lr . *A. ^cosThetaD 
. *1 . /R. ~2*dxE*dyE*dyD*dzD) ; 


*/.•/.  left  face 
clear  V  Vd  Lp 
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Vd(l, : , =  l.xd; 

Vd(2,:,:)  =  yDMax*ones (size (1 . xd) ) ; 

Vd(3 , : , :)  =  l.zd; 

V(l, : , :)  =  Vd(l, : , :)  -  Ve(l); 

V (2 , =  Vd(2 , : , -  Ve(2); 

V (3 , =  Vd(3 , : , -  Ve(3); 

Lp ( 1 , : , : )  =  Lo (1) *ones (size(V(l , :,:))); 

Lp(2,:,:)  =  Lo (2) *ones (size ( V ( 1 ,:,:))) ; 

Lp(3,:,:)  =  Lo (3) *ones (size ( V ( 1 ,:,:))) ; 

R  =  sqrt (V(l ~2  +  V(2,:,:).~2  +  V(3, : , : ) . ~2) ; 
cosThetaE  =  V(3,:,:)./R; 
cosThetaD  =  V(2,:,:)./R; 

A  =  ref lection(nl ,n2 , cosThetaD) ; 

I  =  cookTorranceVL(V,Lp) ; 

Lr  =  A. *1 . *cosThetaE; 

PhiLeft  =  PhiLeft  +  squeeze (Lr . *A . *cosThetaD 
. *1 ./R. ~2*dxE*dyE*dxD*dzD) ; 


*/.•/.  back  face 
clear  V  Vd  Lp 

Vd(l,:,:)  =  -xDMax*ones (size (f .yd)) ; 
Vd(2 , : , : )  =  f .yd; 

Vd(3, : , : )  =  f .zd; 

V(l, : , :)  =  Vd(l, : , :)  -  Ve(l); 

V (2 , :,:)  =  Vd(2, : , :)  -  Ve(2); 

V (3 , :,:)  =  Vd(3, : , :)  -  Ve(3); 

VBack  =  V; 


Lp ( 1 , : , : )  =  Lo (1) *ones (size(V(l , :,:))); 

Lp(2,:,:)  =  Lo (2) *ones (size(V(l ,:,:))) ; 

Lp(3,:,:)  =  Lo (3) *ones (size(V(l ,:,:))) ; 

R  =  sqrt (V(l ,:,:). ~2  +  V(2,:,:).~2  +  V(3, : , : ) . ~2) ; 
cosThetaE  =  V(3,:,:)./R; 
cosThetaD  =  -V(l ,:,:). /R; 

A  =  ref lection(nl ,n2 , cosThetaD) ; 

I  =  cookTorranceVL(V,Lp) ; 

Lr  =  A. *1 . *cosThetaE; 


73 


PhiBack  =  PhiBack  +  squeeze (Lr . *A . *cosThetaD 
. *1 . /R. ~2*dxE*dyE*dyD*dzD) ; 


U  right  face 
clear  V  Vd  Lp 
Vd(l , : , : )  =  l.xd; 

Vd(2,:,:)  =  -yDMax*ones (size (1 . xd) )  ; 
Vd(3 , :  ,  =  l.zd; 

V(l, : , :)  =  Vd(l, : , :)  -  Ve(l); 

V (2 , =  Vd(2 , : , -  Ve(2); 

V (3 , =  Vd(3, : , -  Ve(3); 


Lp ( 1 ,  :  ,  : )  =  Lo (1) *ones (size(V(l , :,:))); 
Lp(2,:,:)  =  Lo (2) *ones (size(V(l ,:,:))) ; 
Lp(3,:,:)  =  Lo (3) *ones (size(V(l ,:,:))) ; 


R  =  sqrt (V(l ~2  +  V(2,:,:).~2  +  V(3, : , : ) . ~2) ; 
cosThetaE  =  V(3,:,:)./R; 
cosThetaD  =  -V(2, : , : ) . /R ; 

A  =  ref lection(nl ,n2 , cosThetaD) ; 

I  =  cookTorranceVL(V,Lp) ; 

Lr  =  A. *1 . *cosThetaE; 

PhiRight  =  PhiRight  +  squeeze (Lr . *A. *cosThetaD 
. *1 . /R. ~2*dxE*dyE*dxD*dzD) ; 


end 


end 


disp(  [’ Completed:  ;  num2str (K/S (1) *100,  ; °/00 . 4g; )  ]) 

end 

PhiVec  =  [sum(sum(PhiTop) ) ;  sum(sum(PhiLeft) ) ; 

sum(sum(PhiFront) ) ;  sum (sum (PhiRight) ) ;  sum(sum(PhiBack) )  ]; 

PhiStruct=  []  ; 

PhiStruct.PT  =  PhiTop; 

PhiStruct.PL  =  PhiLeft; 

PhiStruct.PF  =  PhiFront; 

PhiStruct.PR  =  PhiRight; 

PhiStruct.PB  =  PhiBack; 

PhiStruct . VBack  =  VBack; 

PhiStruct . VFront  =  VFront; 
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m(l)  =  max(max(PhiTop) ) ; 
m(2)  =  max(max(PhiLeft) ) ; 
m(3)  =  max(max(PhiFront) ) ; 
m(4)  =  max(max(PhiRight) ) ; 
m(5)  =  max(max(PhiBack) ) ; 

M  =  1; 

if  0 

figure (1) 
subplot (3, 3, 5) 

image (rot90 (PhiTop, 1)/M*64) ;  title ( ’top’ ) ;  axis  image 
subplot (3, 3, 4) 

image (rot90(PhiLeft ,3)/M*64) ;  title (’left’)  ;  axis  image 
subplot (3, 3, 2) 

image(rot90(PhiFront,4)/M*64) ;  title ( ’ front ’) ;  axis  image 
subplot (3, 3, 6) 

image (rot90 (PhiRight , 1) /M*64) ;  title ( 'right ’) ;  axis  image 
subplot (3, 3, 8) 

image(rot90(PhiBack,2)/M*64) ;  title( ’back’ ) ;  axis  image 
%f igure 

%image (squeeze (IPrime) ) 

end 
if  1 

a  =  size (PhiLeft) ; 
b  =  size (PhiFront) ; 

It  =  [-l*ones (a(l) ,  b(l))  f liplr(PhiFront)/M  -l*ones (a(l) , 
b(l) )]  ; 

Im  =  [rot90 (PhiLeft, 1)/M  fliplr (rot90(PhiTop, 1) )/M 
f liplr(rot90(PhiRight , 1) )/M] ; 

lb  =  [-l*ones (a(l) ,  b(l))  fliplr (f lipud(PhiBack) ) /M  -l*ones (a(l) , 
b(l) )]  ; 

I  =  [It ;  Im;  lb] ; 
figure (2) 

PhiStruct . Image  =  I ; 

imagesc (rot90 (PhiStruct . Image) ) 

axis  image; 

end 


pC  =  sum(sum(PhiTop) )  +  sum(sum(PhiLeft) )  +  sum(sum(PhiFront) ) 
+  sum (sum (PhiRight ) )  +  sum(sum(PhiBack) ) ; 
pE  =  pi*2*xEMax*2*yEMax; 
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eta  =  pC/pE; 

’/figured) 

’/subplot  (3, 3, 9) ; 

[R,DHR]  =  cookTorrance (45, theta) ; 

S  =  eta/DHR 

PhiStruct.xE  =  xe; 

PhiStruct.yE  =  ye; 

end 

end 

function  [A  Rt  Rs  Rp]  =  ref lection (nl ,  n2,  cosTheta) 
sinTheta  =  sqrt(l  -  cosTheta. "2) ; 

Rs  =  abs(  (nl*cosTheta  -  n2*sqrt(l  -  nl/n2*sinTheta) ) 

. / (nl*cosTheta  +  n2*sqrt(l  -  nl/n2*sinTheta) )  ).~2; 
Rt  =  abs(  (nl*sqrt Q-nl/n2*sinTheta)  -  n2*cosTheta) 

. / (nl*sqrt (l-nl/n2*sinTheta)  +  n2*cosTheta)  ).~2; 

Rp  =  (Rs+Rt)/2; 

A  =  1-Rp; 

end 
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